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MEMORANDUM 


Subj : Remote  sensing  of  ocean  significant  wave  height  by  forward 
scattering:  Examples  from  L-band  data 

BACKGROUND 


Forward  scattered  electromagnetic  (EM)  energy  from  the  ocean  is 
generally  detrimental  to  low-angle  tracking  radars,  guidance  systems, 
and  communication  links,  although  positive  use  is  made  of  it  in  some 
height-finding  systems.  In  both  cases,  a thorough  understanding  of  the 
interaction  of  the  EJ1  waves  with  water  waves  is  needed  to  analyze  their 
performance.  Because  Na\'y  applications  often  include  conditions  beyond 
those  of  previous  data,  the  original  purpose  of  tliis  program  was  to 
extend  the  range  of  microwave  frequency  to  1.3  GHz  cL-band)  where  there 
is  a scarcity  of  data.  In  the  course  of  the  data  analysis,  however,  a 
different  application  of  the  data  emerged,  namely,  that  of  remote 
sensing  of  ocean  wave  characteristics. 

A first  series  of  measurements  in  1973  provided  forward-scattered 
data  from  the  slightly  rough  waters  of  Chesapeake  Bay  [1],  A second 
series  in  1974  Involved  reflection  from  the  rougher  waters  off  the 
Atlantic  Coast  using  the  same  airborne  transmitter  flving  at  constant 
altitude  toward  receiving  antennas  mounted  on  a destroyer  (C).  This  is 
the  third  and  final  report  on  the  analysis  of  the  L-band  data  and  the 
remote  sensing  methods  which  have  been  developed. 

FINDINGS 


The  measured  forward  reflection  coefficients  from  the  radar-tracked 

runs  follow  Ament's  [1953]  theoretical  curve  as  a function  of  apparent 

surface  roughness  [ (o  sinii')/Vl  for  the  small  values  (^O.l)  encountered. 

w 

A remote  sensing  methodology  using  forward  scattering  to  measure 

wave  heights  (for  v'alues  of  roughness  below  0.1)  emerges  bv  recasting 

Ament's  expression  into  a series  of  curves  vs  grazing  angle  iij'3 . with 

standard  deviation  of  the  wave  height  (o  ) as  a parameter.  Plotting 

the  measured  reflection  coefficients  on  ¥his  graph  then  determines  wave 

height  (o  ) directly.  The  average  spread  in  the  values  of  o predicted 

bv  this  method  is  -b  0.12  ft  for  values  of  o frc'm  0.4  to  0.1 'ft.  The 
' “ w 

overall  average  error  in  o„  for  six  runs,  compared  to  the  laser  protil- 

ometer  wave  heights,  is  + 5*,  an  Indication  of  the  possible  accviracv 

of  the  method.  An  FM-sweep  method  is  proposed  whic!:  should  improve 

this  accuracy. 

Procedures  have  been  developed  for  possible  remote  sensing  of  the 
two-dimensional  correlation  lengths  of  the  ocean  surface.  These  will 
be  presented  in  a separate  report.  A brief  summarv,  however,  is  given 
in  an  internal  memorandum  [22], 


RESEARCH  IMPLICATIONS 


R€inote  sensing  of  the  ocean  "significant  wave  height"  (Ac  ) by 
forward  reflection  Is  shown  to  be  possible  with  reasonable  accuracy 
for  many  applications,  at  values  of  roughnesses  less  than  0.1. 

Improved  accuracy  should  result  from  increasing  the  frequency  to  10  C.Hz 
and  changing  to  the  FM-sweep  method.  Hie  remotely  sensed  c is 
primarily  that  of  the  peak  spectral  component  of  the  power  Spectral 
curve;  in  the  ocean,  this  is  commonly  the  swell  component. 

There  is  no  established  forward  scattering  theory  for  roughnesses 
above  0.1.  The  implications  are  that  the  two-dimensional  surface 
correlation  lengths  are  involved  at  the  higher  roughnesses  [22].  It  ..s 
necessary  to  resolve  these  basic  questions  experimentally.  Such  an 
experiment  is  discussed  in  the  next  paragraph. 


RECOMMENDATIONS 


.As  there  is  no  adequate  theory  for  the  forward-scattered  coherent 
field  (reflection  coefficient)  for  roughnesses  above  0.1,  an  experiment 
has  been  proposed  to  find  the  basic  mechanisms.  The  different  func- 
tional dependencies  of  shadowing  [20,  IS]  and  spherical  wave  theories 
[17,  5,  15]  and  hypotheses  involving  the  two-dimensional  correlation 
lengths  [22]  would  be  exploited  in  establishing  the  mech.inisms. 

Measurements  would  be  made  under  controlled  two-dimensional  water-wave 
conditions,  with  adjustable  wave  directions,  heights,  lengths,  and 
two-dimensional  surface  correlation  functions.  To  satisfy  these  two- 
dimensional  conditions,  a wave  basin  is  required  (e.g.  , David  Taylor 
Model  Basin).  To  achieve  roughnesses  above  0.1,  a short  EM  wavelength 
of  3.3  mm  (90  GHz)  is  required.  This  frequency  will  also  be  advantageous 
in  extending  the  wavelength  region  of  forward  scattering  data. 
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Remote  Sensing  of  Ocean  Significant  Wave  Height 
by  Forward  Scattering:  Examples  from  L-Band  Data 

1 . INTRODUCTION 


The  measurements  described  previously  [1,  2]  and  In  this  report  were 
originally  carried  out  to  obtain  data  on  the  forward  reflection  coeffi- 
cient of  the  ocean  at  an  electromagnetic  wavelength  of  23,1  cm  (L-band). 
Years  ago,  measurements  of  forward  scattering  of  5.7-  to  35-GHr  micro- 
wave  radiation  from  the  surface  of  the  ocean  led  to  the  establishment  of 
a scattering  model  [3,  4,  5]  for  both  the  coherent  and  Incoherent 
scattered  fields  as  a function  of  "apparent  surface  roughness."* 

Although  this  model  exists.  Navy  applications  often  Include  conditions 
beyond  those  of  the  data  forming  the  model.  For  example,  valid  L-band 
data  (with  ocean  wave  measurements)  at  low  grazing  angles  are  scarce 
[6].  One  purpose  of  this  program  was  to  extend  the  range  of  parameters 
to  a lower  frequency  (1.3  GHz)  so  that  the  model  could  serve  a wider 
range  of  applications. 

Such  data  provide  information  for  analysis  of  a radar  being  used  for 
height  finding  of  aircraft,  for  design  of  low-angle  tracking  radars,  for 
use  in  simulators  of  low-angle  tracking  of  aircraft  and  cruise  missiles, 
for  analysis  of  over-ocean  communication  circuits,  IFF,  and  missile 
guidance,  and,  in  general,  for  any  Navy  application  in  which  an  antenna 
beam  is  directed  over  water  such  that  some  water-scattered  energy  is 
received. 

A new  application  emerged  from  the  analysis  of  the  data,  namely, 
remote  sensing  of  ocean  wave  parameters.  Determination  of  the  "signifi- 
cant wave  height"  is  demonstrated  in  this  report  with  the  L-band  data. 
The  measured  L-band  forward  reflection  coeff iclents,  when  plotted  on  a 
graph,  yield  the  "significant  wave  height"  (4o  ) directly.  In  addition, 
procedures  are  being  b*orked  out  for  possible  remote  sensing  of  the  two- 
dimensional  correlation  lengths  of  the  ocean  surface.  These  will  be 
presented  in  a future  report. 

A first  series  of  measurements  in  1973  provided  data  of  reflection 
from  the  slightly  rough  waters  of  Chesapeake  Bay  as  well  as  a system 


* Apparent  surface  roughness  is  defined  here  as  the  "phase  roughness"  of 
the  surface,  i.e.,  (o  3in<t>)/\,  where  o is  the  standard  deviation  of  the 
surface  height  variatYons,  C'  is  the  grazing  angle,  and  \ is  the  electro- 
magnetic wavelength.  This  quantity  can  be  converted  to  phase  in  radians 

by  multiplying  by  2*:. 

Note:  Manuii'ript  tubmitted  JMiuarv  U1.  U179. 
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test  []].  A fixed  receiver  on  a cliff  overlooking  the  Bay  received 
energy  from  a transmitter  in  an  aircraft  flying  at  constant  altitudes. 

A second  series  involved  reflection  from  the  rougher  waters  off  the 
Atlantic  Coast  using  the  same  airborne  transmitter  but  with  the 
receiving  antennas  mounted  on  a stabilired  platform  of  a destroyer.  Six 
days  of  data  were  obtained  during  the  period  February  to  .■Vprll  1*^74,  and 
a preliminary  report  was  written  [2].  This  is  the  final  report  on  these 
data. 

The  report  covers  two  principal  topics,  u)  remote  sensing  of  ocean 
wave  parameters,  and  (b)  detailed  analysis  of  the  L-band  data.  Those 
readers  who  may  be  Interested  only  in  the  remote  sensing  aspects  should 
see  Sections  S.2,  6.1,  6.2,  6.3,  and  6.6.  In  reading  Section  5.2, 
concentrate  principally  on  the  results  from  radar-tracked  Missions  17 
and  18;  in  Missions  13,  14,  and  15  the  ship's  radar  was  not  working,  and 
optical  tracking  presumably  produced  wild  swings  in  the  data. 

2.  DESCRIPTION  OF  TESTS  .ANP  EQUIPMENT 

The  measurement  system  was  relatively  simple  and  consisted  of  a CW, 
airborne  transmitter  and  a shipboard  receiver.  Both  the  direct  ray 
from  the  transmitter  and  the  surface-reflected  signal  entered  the  wide- 
beam  receiving  horns.  .As  the  aircraft  approached  the  ship,  the  received 
signal  traced  out  the  usual  interference  (or  lobe)  pattern  formed  by  the 
direct  and  reflected  signals,  which  combine  at  a continuously  increasing 
phase  angle  corresponding  to  the  increasing  path  length  difference 
between  direct  and  specularly  reflected  rays.  .As  details  of  the  equip- 
ment and  experiment  are  given  in  (2l,  only  a brief  summary  is  Included 
here  for  continuity. 

The  L-band  Cm'  transmitter  was  mounted  within  an  oval  shaped  pod 
suspended  from  the  outboard  port  bomb  rack  of  the  NRL  S-2D  aircraft. 

The  antenna,  a 15.5-dB  standard  gain  rectangular  hern,  liad  its  central 
axis  depressed  8*  below  the  horizontal  to  Improve  ilUimination  of  the  >.] 

water.  .Although  preflight  selection  of  either  horizontal  or  vertical 
polarization  was  possible,  only  horizontal  polarization  was  used  in 
these  measurements.  .Antenna  beamwidth  was  approximately  32*  in  both  ii 

planes  with  patterns  shown  in  [2]. 

Two  receiving  horns  (15.5  dB,  32*  beamwidth^  one  horizontally 
polarized  and  one  vertically  polarized,  were  mounted  50  feet  above  the 

water  on  the  antenna  yoke  of  the  Mark  25  fire  control  radar'gun  director  ^ 

on  each  of  tw  destroyers,  the  I’SS  JOHNSTON  (DD-821)  and  the  I’SS  LAFFEY 
(DD-''24).  This  mounting  point  was  cross-level  stabilized  to  minimize 
the  effect  of  roll  and  pitch  on  the  received  signal. 

Ship  motion  was  recorded  aboard  the  I'SS  JOHNSTON.  The  outputs  of 
the  ship's  gyro  system  were  digitally  recorded  on  a separate  time- 
synchronized  tape.  The  ship  motion  data  were  needed  for  receiver 


antenna  orientation  corrections  for  interpretation  of  the  Interference 
patterns.  Wind  speed  and  direction  were  obtained  aboard  ship  from  an 
anemometer  ^13  feet  above  sea  level. 

Data  were  taken  as -the  S-2D  aircraft  flew  directly  toward  the  ship. 
Either  the  Mark  25  radar  or  an  optical  tracker  cracked  Che  incoming 
aircraft  to  point  the  receiving  horns  and  to  provide  range  data.  Only 
two  flight  altitudes  were  used  (1000  and  2500  feet)  to  enable  repeat 
runs  at  the  same  altitude  in  an  attempt  Co  obtain  a better  statistical 
sample.  The  pilot  estimated  that  he  maintained  the  altitude  to  within 
+ 20  feet  of  these  values. 

Ocean  wave  heights  were  obtained  by  an  airborne  laser  profllomeCer 
(7].  For  the  series  of  tests  with  the  USS  JOHNSTON,  a separate  Piper 
Navaho  aircraft  (provided  by  NOVA  University)  carried  the  laser.  In 
the  later  test  with  USS  LAFFEY,  the  laser  profilometer  was  mounted  in 
the  same  NRL  S-2D  aircraft  which  carried  the  transmitter. 

Aerial  photographs  were  also  taken  at  an  altitude  of  1000  feet  by 
a camera  in  Che  transmitter  pod  for  analysis  by  the  Stilwell  method 
[8,  9].  The  wide-angle  photographs  needed  for  density  normalization  in 
this  analysis  were  taken  from  shipboard. 

3.  EXPERIMENTAL  DATA 


3.1  Summary  Tables 

Data  were  collected  at  different  locations  off  the  Atlantic 
Coast  from  Florida  to  Virginia  on  six  different  days  during  February  to 
April  1974,  as  listed  in  Table  I.  All  transmitter  aircraft  altitudes 
were  1000  feet,  excpet  for  those  runs  marked  by  an  asterisk.  The 
headings  of  the  transmitter  aircraft  relative  to  the  wind  direction  are 
given  for  each  run,  where  known,  in  Table  II.  The  transmitting  horn 
was  oriented  for  horizontal  polarization  in  level  flight  for  these  six 
days  of  data.  All  data  in  this  report  are  for  horizontal  polarization, 
except  for  Fig.  12  from  an  earlier  report  [1]. 

In  the  shipboard  receiving  system,  a coaxial  diode  switch 
alternately  connected  the  single  channel  receiver  input  to  the  two 
orthogonally  polarized  (horizontal  and  vertical)  receiving  horns.  The 
switch  was  driven  at  a 100-Hz  rate.  This  situation  is  denoted  by  the 
symbol  HV  in  Table  I.  In  the  remaining  runs  (denoted  by  H)  only  hori- 
zontal polarization  was  recorded.  The  data  were  digitally  recorded  on 
magnetic  tape  at  a sampling  rate  of  100  Hz,  synchronized  to  the  r-f 
switching.  The  times  of  the  data  runs  are  given  in  Appendix  A. 

3.2  Examples  of  Interference  Pattern  Data 

Analog  strip  chart  recordings  were  also  made  to  monitor  the 
signals  in  real  time.  Figure  1 is  an  example  of  the  analog  charts 
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showing  d continuous  Interfersnc*  pattsrn  for  an  aircraf c-cransmitter 
alclcuds  of  1000  feat.  A dB-scala  of  relative  signal  strength  has  been 
placed  on  the  ordinate.  On  the  abscissa,  elapsed  time  increases  to  the 
right  as  the  aircraft  carries  the  transmitter  toward  the  receiver.  A 
time  scale  and  approximate  ranges  in  kiloyards  are  marked  on  the  chart. 

Ac  ranges  less  chan  three  kiloyards,  Che  record  is  unusable  because  the 
specular  ray  is  "sliding  down  the  side"  of  Che  antenna  beam. 

In  comparison,  the  regularity  of  an  interforence  pattern 
measured  over  the  calmer  waters  of  Chesapeake  Bay  [1,  Fig.  7]  is  striking 
compared  to  the  over-ocean  results  of  Fig.  1. 

Instead  of  the  condensed  analog  trace  of  Fig.  1,  sections  of 
expanded,  callbratad-computer  plots  are  given  in  Fig.  2.  These  are 
different  sections  of  one  run,  i.e.,  during  one  incoming  aircraft  pass. 

The  signal-strength  scale  is  in  dB  below  an  arbitrary  level.  The 
elapsed-cime  scale  is  one  second  for  one  inch  on  the  original  computer 
plot  for  the  several  minute  runs.  All  the  small  spikes  and  the  sporadic 
large  spikes  in  Fig.  2 are  noise  and  should  be  ignored.  Various  t\T<s 
of  Interference  minima  are  encountered,  as  illustrated  bv  the  samples 
in  Fig.  2.  These  Include  the  sharp,  clean,  deep  type  in  Fig.  2a,  the 
broad  type  in  Fig.  2b,  .and  smaller  apparent  minima  on  the  slopes  of  larger 
minima  in  Figs.  2c  and  2d.  These  are  discussed  in  Section  4. 

4.  DATA  REDUCTION 


4. 1 Calculated  Lobe  positions 

The  calculated  positions  of  the  interference  pattern  minima 
are  used  to  select  the  proper  minima  and  reject  spurious  minima  which 
arise  from  causes  to  be  discussed.  The  simplest  procedure  will  be 
used  to  calculate  the  positions  of  the  maxima  and  minima  by  assuming 
a,fla^,  smooth  surface,  no  atmosphere  (or,  a homogeneous  one)  and  that 
Z‘.  Z‘  <<  R“,  where  Z.  and  are  the  receiver  and  transmitter  heights, 
respectively , and  R is  the^range.  Then^tlje  difference  between  the 
direct  ray  p.ath  length,  ^ [R“ ,+  and  the  specularly  reflected 

ray  path,  [R*  (Z.+Z^^) * 1 . is  approximately : 

AL  a 2Z^Z,/R. 


Maxima  and  minima  occur  when  AL  ■ n\,'2,  where  n is  an  integer.  Then 
the  locations  of  maxima  and  minima  are  given  approximately  bv: 


Table  III  lists  calculated  values  of  R for  various  n to  compare  with 
Fig.  1.  The  values  are  alternat'ly  maxima  and  minima,  but  without 
designation.  Normally,  maxima  would  occur  for  odd  values  of  n because 
the  reflection  coefficient  for  horizontal  polarization  is  near  -1.  As 
these  simplified  calculations  do  not  include  an  atmosphere,  antenna 
differential  phase  shifts,  etc.,  the  calculated  values  can  be  shifted  as 
a group  to  obtain  the  best  match. 

Refinements  to  these  simple  calculations  could  incorporate  a 
standard  atmosphere,  (4/3  radius  earth  approximation),  and  spherical 
earth  geometry,  if  warranted.  Programs  for  these  calculations  have  been 
worked  out  by  Blake  [10], 

Two  of  the  best  fits  between  calculated  and  actual  minima 
positions  are  for  Run  2 of  Mission  14,  given  in  Table  IV,  and  Run  5 of 
Mission  15  in  Table  V.  Last  in  Table  VI  is  an  example  of  an  inferior 
fit.  The  best  match  to  the  even-integer  calculated  ranges  (the  ones 
matching  the  other  tiro  runs  so  well)  is  obtained  when  n “ 6 at  the  most 
distant  observed  minimum  of  15.34  kyds.  The  match  to  the  odd  set  of 
integers,  beginning  with  n - 7,  is  just  as  good,  but  with  the  added 
point  that  the  differences  between  calculated  and  observed  positions 
are  constant  at  - 2 kyds  (after  the  first  minimum).  How  such  a differ- 
ence could  have  arisen  is  not  clear.  Could  a bias  have  existed  in  the 
optical  tracker,  in  reading  ranges  from  the  TACAN  system,  etc.,  during 
the  first  run  of  the  day?  The  correspondences  to  the  sets  beginning 
with  n “ 5 or  n ■ 8 are  worse  than  those  shown  for  n ■ 6 and  n " 7. 
kTien  the  fit  of  the  absolute  positions  is  poor,  as  in  Table  VI,  then 
first  differences  are  used  to  help  select  the  proper  minima  and  reject 
the  extraneous  minima  (the  latter  are  discussed  below). 

Radar  tracking  ranges  became  available  for  reducing  the  data 
when  the  Mk  25  radar  first  became  operational  in  Mission  16.  Problems 
are  encountered,  however,  with  calculated  lobe  positions  and  with  the 
specular  ray  "sliding  off  the  beam"  at  an  unlikely  value  of  the  grazing 
angle.  It  is  found  that  the  radar  ranges  in  Missions  16  and  17  are 
greater  chan  the  optical  tracker  ranges  by  10,000  ft  and  11,000  ft, 
respectively.  Figure  3 shows  this  discrepancy  for  one  run  of  Mission 
17.  In  Mission  18  on  a different  ship,  and  hence  with  a different 
radar,  the  radar  and  optical  ranges  agree  (Fig.  4).  The  trouble  could 
have  been  either  in  the  radar  (which  had  to  be  repaired)  or  in  the 
range  conversion  unit  which  converted  time  delay  into  digital  range.  In 
Mission  18,  the  fact  that  the  radar  and  optical  ranges  agree  tends  to 
exonerate  the  range  conversion  unit.  Hence,  the  interference  pattern 
data  of  Missions  13  through  17  are  computed  using  optical  tracker  ranges, 
while  radar  ranges  are  used  in  Mission  18,  as  they  are  recorded  directly 
on  cape. 

Some  possible  causes  for  disagreement  between  the  computed  and 
actual  patterns,  are  as  follows. 
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4.1,1  Incoherent  Field 


The  incoherent  field  (of  random  magnitude  and  phase)  is 
produced  by  scattering  of  the  electromagnetic  waves  at  the  rough  sea 
surface.  This  incoherent  field  adds  to  all  parts  of  the  Interference 
pattern,  but  produces  the  greatest  visible  effect  at  the  pattern  minima, 
where  the  signal  is  weak,  being  the  difference  of  the  direct  and  coherent 
fields.  Thus,  although  the  airborne  transmitter  in  constant  altitude 
flight  has  to  "fly  through"  each  interference  minimum,  there  is  the 
possibility  that  a correct  minimum  may  not  be  the  deepest  of  the  nearby 
spurious  (or  real)  minima.  The  reason  is  that  the  incoherent  field  may 
be  additive  at  the  particular  instant  when  the  aircraft  "flies  through" 
it,  thus  raising  the  depth  of  the  minimum.  Background  information  on 
the  incoherent  field  is  available  [3,  4,  5],  and  some  preliminary  results 
for  the  Chesapeake  Bay  data  have  been  given  [1]. 

4.1.2  Vertical  Motion  of  Aircraft 

The  calculated  lobe  patterns  assume  that  the  aircraft  trans- 
mitter flies  at  a constant  altitude.  If  it  does  not,  the  positions  of 
the  lobes  will  be  altered.  Table  VII  tabulates  values  of  height  changes, 
AZ2,  which  will  produce  phase  changes  of  180®,  90°,  45®,  and  30®.  At 
the  longer  range  of  75,000  feet,  all  values  of  AZ2  are  too  large  to  be 
of  concern.  At  the  shorter  range  of  10,000  feet,‘^however , the  required 
values  of  AZ^  for  phase  changes  of  90®  or  less  are  within  the  + 25-foot 
altitude  tolerance  of  the  aircraft.  As  an  example,  if  the  aircraft 
passes  through  a minimum  and  is  45®  in  phase  away  from  It,  a drop,  AZ, , 
of  9 ft  would  return  the  signal  to  the  minimum;  the  recorded  signal  would 
show  an  extra  minimum  riding  on  the  shoulder  of  a minimum. 


Rapid  height  changes  are  not  necessary  for  all  the  effects. 
If  the  aircraft  descends  at  the  proper  constant  angle  (a  multiple  of 
X/2Zi,  as  seen  from  a point  on  the  water  under  the  receiving  horn),  it 
can  "ride  down"  a minimum.  This  behavior  might  be  the  cause  of  the 
broad  minimum  in  Fig.  2b. 

Also,  a AZ^  of  a size  to  shift  from  a maximum  to  a minimum 
is  not  necessary  to  cause  good-sized  minima.  Table  VIII  lists  signal 
level  changes  for  various  phase  angles  from  a minimum.  For  a 30°-phase 
change  from  the  minimum,  the  signal  will  change  14  dB;  this  is  for  a 
AZ-  of  only  6 ft  at  the  shortest  range  of  10,000  ft.  Instead  of  37  feet 
for  the  full  180®  change.  In  Fig.  2c,  from  54  to  56  seconds,  the  air- 
craft may  be  bobbing  up  and  dovm  one  or  two  feet,  possibly  from  turbu- 
lence, to  produce  the  signal  fluctuations.  A drop  of  six  feet  is  then 


sufficient  to  reach  the  minimum  of  58  seconds. 

In  Fig.  2d,  there  are  three  approximately  constant-signal 
regions,  below  the  signal  levels  of  the  maxima,  where  the  aircraft  may 
be  flying  slightly  downward  at  the  proper  angle  (as  described  for  Fig. 
2b)  to  remain  at  a constant  signal  level  position  in  the  interference 


pattern.*  (These  three  regions  are  from  6 to  7 secs,  8 to  8 1/2  secs, 
and  9 to  10  secs.)  A six-foot  drop  in  altitude  would  then  account  for 
the  minimum  at  11  secs. 

4.1.3  Vertical  Motion  of  First  Fresnel  Zone 

Even  if  the  aircraft  flies  at  a constant  altitude,  the  swell 
moves  sizeable  areas  of  water  up  and  down.  Usually  the  Fresnel  zones 
are  long  enough  to  extend  over  several  swell  wavelengths,  but  at  the 
shortest  ranges  this  may  no  longer  be  true.  Estimates  of  this  effect 
can  be  made  by  the  following  calculations: 

22,  Z, 

Again,  AL/\  » 


Tr  V'l'- 


For  a water  level  change,  dZ  « dZ,  - dZ*, 

1 J dZ 


As  an  example,  for  ■ 50  ft,  ■ 1000  ft,  R 


20,000  ft 


f ■ 1.3  GHZ,  — ^ * 0.13*5  * 50®  phase  change  per  foot  of  dZ. 

For  2^  " 2500  ft,  * 0.336  ^ 120®  ph.ise  change  per  foot 

of  dZ. 


Thus,  if  the  swell  could  lift  the  entire  1st  Fresnel  zone  up  and  down 
only  one  foot,  the  signal  would  move  through  a good  fraction  of  the 
Interference  minimum.  Some  example  dimensions  of  the  1st  and  10th  Fresnel 
zones,  calculated  by  the  equations  in  Kerr  [11],  are  given  In  Table  IX 
for  aircraft  altitudes  of  1000  and  2500  ft.  Where  \ Is  the  distance  of 
the  center  of  the  1st  Fresnel  zone  from  the  receiver*,  i is  the  length  of 
the  semi-major  axis  of  the  ellipse  (along  the  line  between  receiver  and 
transmitter),  and  y is  the  semi-minor  axis.  Thus  for  2^  “ 2500  ft  and 
R ■ 20,000  ft,  the  1st  Fresnel  zone  is  centered  at  422  ft  (from  the 
receiver)  and  extends  + 137  ft  from  this  point,  or  from  285  to  559  ft. 


For  deep-water  waves. 


\ 

w 


where  \ is  the  wavelength  In  feet  if  g Is  the  acceleration  of  gravity 

in  feetysecond^,  .and  T is  the  period  In  seconds.  Ex.ample:  T - 7 secs., 

\ ■ 250  ft,  T ■ 8 secs.,  \ ■ 328  ft,  and  T ■ 9 secs.,  \ ■ 415  ft. 
w w w 

* Removing  the  1/R  factor  from  the  discussion. 
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If  the  transmitter  aircraft  files  parallel  to  the  wave  crests, 
the  narrow  1st  Fresnel  zone  fits  between  the  crests  and  could  rise  and 
fall  as  an  entity;  short-crested  waves  would  break  this  up  to  some  extent. 

If  the  aircraft  files  either  up-  or  downwind  courses,  only 
the  smallest  Fresnel  zones  in  Table  IX  might  qualify.  For  example,  at 
Z2  ■ 2500  ft  and  R ■ 20,000  ft,  the  1st  Fresnel  zone  of  274-foot  length 
barely  fits  between  the  crests  of  the  S-second  wave,  and  obviously  a 
crest  could  lift  only  part  of  this  zone.  The  9-second  wave  of  415  ft 
wavelength, lifting  and  dropping  a significant  portion  of  the  1st  Fresnel 
zone, could  cause  at  least  a partial,  or  small  Interference  minimum.  It 
appears  that  this  phenomenon  would  only  be  operative  at  the  shortest 
ranges  (<  20,000  ft)  for  an  aircraft  height  of  2500  ft  when  the  course 
Is  up-  or  downwind.  The  lifting  and  lowering  of  the  Fresnel  zone  would 
occur  at  the  slow  swell  period,  and  thus  could  account  for  several 
minima. 


4.1.4  Shipboard  Antenna  Height  Changes 

The  receiving  horns,  mounted  50  ft  above  the  water,  are  sub- 
ject to  changes  In  height  as  the  ship  rolls  and  pitches.  The  effect  of 
this  height  change,  Is  calculated  as  In  Section  4.1.2  for  vertical 

motion  of  the  aircraft. 

r-  S-i-R 

‘^‘'1  2Z^  n 

Let  n ■ 2,  or  a 180*  phase  change,  and  R ■ 10,000  ft;  If  Z^  ■ 2500  ft 
then  AZ^  ■ 0.75  ft;  If  • 1000  ft,  then  AZj^  “ 1.9  ft. 

Thus,  for  an  altitude  of  2500  ft  at  the  shortest  ranges,  the 
required  change  In  AZ^  Is  easily  achieved  with  ship  rolls  of  less  than 
10*.  The  maximum  change  of  height  occurs  four  times  during  a roll 
period,  or  with  typical  destroyer  roll  periods  of  9 seconds,  every 
2 1/4  seconds.  Hence,  this  mechanism  could  be  a likely  cause  of  several 
additional  minima  at  the  shorter  ranges.  Roll  and  pitch  were  recorded 
on  the  USS  JOHNSTON,  but  these  records  have  not  been  reduced. 

4.2  Plots  of  ’’Instantaneous  Reflection  Coefficient"  vs  Grazing  Angle 

If  the  reflecting  surface  were  a perfectly  smooth  plane,  the 
Interference  pattern  would  be  smooth,  because  there  would  be  no  fluctu- 
ating Incoherent  scattered  field.  Signal  maxima  occur  when  the  direct 
and  coherent  reflected  ravs  add  in  phase,  and  minima  occur  when  they 
are  130*  out  of  phase. 

The  value  of  the  amplitude  reflection  coefficient,  r,  is  then 
calculated  from  adjacent  maxima  and  minima  by  the  usual  expression 
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where  E^and  E are  the  electric  field  strengths  at  the  maxima  and  minima, 
respectively.  At  longer  ranges  where  the  spherical  earth  becomes  impor- 
tant, r has  to  be  multiplied  by  a divergence  factor  [11,  12].  These 
spherical  earth  effects  are  not  important  at  the  ranges  and  altitudes 
used  in  these  tests.  When  the  surface  is  not  smooth,  this  simple 
expression  is  no  longer  correct  because  of  the  addition  of  the  random 
Incoherent  field;  a statistical  treatment  then  has  to  be  used.  Never- 
theless, eq.  (1)  is  applied  to  the  successive  maxima  and  minima  of 
several  runs  to  form  what  might  be  called  an  "instantaneous  reflection 
coefficient,"  The  purpose  is  to  show  the  scatter  of  the  data  points — 
the  results,  however,  are  startling.  Instead  of  a random  scatter  of 
points  about  a gradually  decreasing  curve  of  "instantaneous  reflection 
coefficient"  as  the  grazing  angle  Increases,  the  data  points  depict  a 
series  of  oscillations,  as  shown  in  Fig.  5 for  Run  9 of  Mission  13.  The 
two  large  circle  points  are  for  smoothing  over  1®  grazing  angle  intervals, 
as  described  later.  Table  X lists  the  data  points  from  which  Fig.  5 is 
plotted. 

The  theoretical  curve  of  the  coherent  field  reflection 
coefficient  versus  the  apparent  surface  roughness  parameter  (a^siniJ))/X 
is  a Gaussian  shape,  as  derived  by  Ament  [13].  (See  plot  later  in 
Fig.  13.)  Assuming  that  a is  constant  during  one  data  run  (an  Inbound 
aircraft  flight),  and  wlth'^X  fixed,  the  abscissa  variable  becomes  sini/, 
or  Just  ^ for  small  angles  as  in  Fig.  5. 

The  angular  separation  between  the  four  peaks  of  "instantaneous 
coefficient"  in  Fig.  5 is  approximately  a constant  of  0.40*  + 0.05*.  The 
time  Intervals  between  the  successive  peaks  continually  decrease  as  the 
aircraft  comes  in  (Table  X). 

A plot  for  another  run  (Run  12)  on  Mission  13  shows  similar 
behavior  (Fig.  6).  The  angular  separation  between  peaks  of  the 
"Instantaneous  reflection  coefficient"  is  again  approximately  a constant, 
being  ~0.6*  + 0.1*,  and  the  time  intervals  between  peaks  also  decrease 
as  the  grazing  angle  increases.  Table  XI  contains  the  data  for  Fig.  6. 

The  plot  for  Run  8 of  Mission  15  shows  a remarkable  series  of 
six  peaks,  the  last  five  of  which  have  a regular  spacing  of  0.87*  + 0.08* 
and  decreasing  spacing  in  time  (Fig.  7 and  Table  XII),  as  do  the  two 
runs  of  Mission  13. 

On  these  two  days  (Missions  13  and  15),  only  manual  optical 
tracking  of  the  aircraft  was  used  because  the  Mk  25  automatic  tracking 
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radar  had  a fault  and  was  inoperative.  The  suggestion  was  made*  that 
the  manual  optical  tracking  might  be  responsible  for  the  sequential 
oscillations  in  the  "instantaneous  reflection  coefficient"  plots.  The 
assumed  mechanism  would  be  that  as  the  receiver  antenna  pointing  over- 
shot and  pointed  above  the  aircraft,  the  reflected  specular  "ray"  would 
slide  down  the  side  of  the  antenna  beam,  resulting  in  too  small  an 
"Instantaneous  reflection  coefficient."  Conversely,  as  the  receiver 
antenna  pointing  overshot  in  the  other  direction  and  pointed  below  the 
aircraft,  the  direct  ray  would  slide  down  the  side  of  the  beam  to  cause 
too  large  an "instantaneous  reflection  coefficient." 

Figures  8 and  9 (and  Tables  XIII  and  XIV)  made  to  test  this 
hypothesis  show  one  run  on  each  of  the  last  two  days  (Missions  17  and 
18)  when  automatic  radar  tracking  of  the  aircraft  is  used  to  point  the 
receiving  horns.  Indeed,  Fig.  8 looks  like  a marked  departure  from 
Figs.  5,  6,  and  7 in  that  the  oscillations  are  "suppressed"  and  do  not 
begin  until  5*.  Figure  9 (Run  3 of  Mission  18)  also  looks  similar  to 
Fig.  8 in  that  the  oscllllations  are  "suppressed"  until  3 3/4®.  Thus, 
surprisingly,  both  Figs.  8 and  9 with  automatic  radar  tracking  show  some 
oscillatory  behavior  (albeit  at  higher  grazing  angles,  or  shorter  ranges) 
similar  to  the  manual  optical  tracking  data  in  Figs.  5 to  7 . As  the 
oscillations  are  apparently  reduced  at  the  longer  aircraft  ranges  (lower 
grazing  angles)  by  the  radar  tracking,  it  could  be  inferred  that  perhaps 
the  tracking  loop  or  stabilization  system  is  not  adequate.  For  example, 
is  it  possible,  as  the  angular  tracking  rates  Increase  at  the  shorter 
aircraft  ranges,  that  the  ship's  roll  and  pitch  are  not  adequately 
compensated?  If  so,  this  would  be  true  for  both  the  two  different 
destroyers  used  in  Missions  17  and  18.  Ship  motion  data  were  recorded 
for  Missions  13  and  17,  but  not  for  Mission  18.  There  has  not  been 
sufficient  time  to  reduce  these  records. 

\^at  can  be  done,  however,  is  to  compare  results  of  aircraft 
data  runs  made  over  Chesapeake  Bay,  with  the  receiver  horn  fixed  on  a 
building  atop  the  cliff  at  the  edge  of  the  Bay.  Thus,  there  will  be  no 
ship  motion  Involved,  but  only  the  motion  of  the  aircraft.  The  next 
three  Figs.,  10,  11,  and  12,  are  data  from  the  Chesapeake  Bay  report  [ll 
with  the  addition  of  connecting  lines  between  the  points  for  easier 
comparison  to  Figs.  5 to  9.  Uliether  any  regularity  still  remains  in 
Figs.  10  and  11  is  doubtful;  there  does  not  appear  to  be  any  regularity 
in  Fig.  12.  In  any  case,  the  elimination  of  the  ship  with  its  tracking 
system  as  a receiving  terminal  has  at  least  greatly  reduced,  if  not 
eliminated,  any  apparent  regularity. 

Whatever  the  cause,  the  large  "instantaneous  reflection 
coefficient"  swings  and  their  regularity  cast  doubt  on  the  validity  of 
the  data,  unless  properly  handled.  In  this  vein,  it  is  reassuring  to 
see  how  well  the  points  using  1®  grazing  angle  interval  smoothing 
(described  later)  appear  visually  to  describe  the  average  of  the  swings 
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in  Figs.  5 to  9.  (Grazing  angle  smoothing  intervals  of  1/2*  are  found 
to  be  Inadequate;  Figs.  5 to  7 show  one  reason  why.)  Thus,  data  results 
which  use  smoothing  intervals,  such  as  the  coherent  field  reflection 
coefficient,  can  be  obtained  from  Missions  13  to  18  with  shipboard 
receiving  terminals;  a meaningful  incoherent  field  cannot  be  obtained, 
however,  because  its  computation  uses  Instantaneous  signal  values,  which 
may  be  contaminated  by  ship  motion  effects. 

4 . 3 Calculation  of  Coherent  Field  Reflection  Coefficient 

4.3.1  Signal  Model 

To  calculate  the  reflection  coefficient  of  the  coherent 
field  from  total  field  data,  it  is  necessary  to  extract  the  incoherent 
scattered  field  from  the  total  field.  Use  can  be  made  of  a statistical 
forward  scattering  model  based  upon  experimental  data  [3,  4]. 

A brief  recapitulation  is  included  here  for  completeness. 

Let  the  instantaneous  total  received  field  be  E^.  The  temporal  envelope 
of  this  field,  fluctuating  at  the  doppler  frequencies,  is  what  is 
detected  and  measured;  for  simplicity  of  notation,  it  will  be  denoted  as 
T.  This  total  field  is  composed  of  a direct  field  D,  a coherent  reflected 
field  (of  magnitude  C and  phase  0,  and  a random  incoherent  field  of 
instantaneous  magnitude  I.  The  strength  of  the  direct  field  varies 
inversely  with  slant  range  from  receiver  to  aircraft  transmitter.  If 
both  the  transmitter  and  receiver  were  fixed  Instead  of  moving,  then  the 
magnitude  C and  phase  0 would  be  constant  in  a time-average  sense.  The 
Incoherent  field  I is  usually  assumed  to  arise  from  a large  number  of 
independent,  random,  surface  scatterers  and  thus  to  be  Rayleigh 
distributed  with  equiprobable  phase.  For  deviations  from  the  Rayleigh 
distribution,  see  [5].  A combined  "signal"  phasor  S is  introduced  as 
S ■ [d  + c|,  the  phasor  sum  of  D and  C.  The  "signal  plus  noise"  statis- 
tics of  S.O.  Rice  [14]  are  used  to  analyze  the^rec^lved  total  field 
(for  a Rayleigh  Incoherent  field),  using  S - [d  + c|  as  the  "signal" 
and  the  Incoherent  field  I as  the  "ncise." 

4.3.2  Effects  of  the  Changing  Aircraft  Range 

In  the  experiments  establishing  the  model,  both  the 
transmitter  and  receiver  positions  were  fixed  during  a data  run,  and  a 
recording  of  the  received  field  was  made  over  a time  interval  (e.g., 

10  minutes)  to  obtain  the  fluctuating  signal  time  average  and  variance. 

In  the  measurements  of  this  report,  however,  the  transmitter  position 

moves  continuously,  consequently  rotating  the  phase  angle  (0)  between  i 

C and  D continuously,  and  giving  rise  to  the  maxima  and  minima  of  the 

signal.  There  are  five  factors  (at  least)  affecting  the  amplitude  of 

the  field  as  the  range  changes.  First,  the  magnitude  of  C changes  in 

a functional  relationship  w;^th  apparent  ocean  roughness,  ^ = (c  sinC')/\, 

as  C(rough)  ■ exp  - 2(2ir*)-,  as  derived  by  Ament  [13].  As  the'^alrborne 

transmitter  approaches  at  constant  altitude,  the  grazing  angle  (if) 
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increases,  and  C (rough)  decreases  according  to  Ament’s  expression. 
Second,  the  smooth  surface  reflection  coefficient,  C (smooth),  varies 
with  grazing  angle  (and  polarization) , as  given  by  the  curves  in 
Figs.  5,3  to  5.6  of  Kerr  [11].  (Other  influencing  factors  such  as 
salinity,  temperature,  surface  contamination,  etc.,  have  to  be  assumed 
to  be  constant  over  the  propagation  path.)  Third,  the  spherical  earth 
divergence  factor  will  be  neglected. 

A fourth  factor  arises  because  of  differences  in  illu- 
mination of  the  surface.  Ament's  [13]  expression  was  derived  assuming 
plane  waves,  whereas  spherical  waves  modified  by  various  antenna  beam 
patterns  are  used  experimentally.  This  effect  cannot  be  expressed  by  a 
simple  multiplication  of  an  antenna  beam  illumination  factor  times 
C(rough),  as  has  been  done  in  the  past.  For  example,  a narrow  beam, 
pointed  horizontally,  will  Illuminate  the  Fresnel  zones  on  the  distant 
side  of  the  specular  point  almost  uniformly,  whereas  the  Fresnel  zones 
on  the  near  side  of  the  specular  point  can  be  illuminated  very  non- 
unlformly.  This  Illumination  changes  as  the  aircraft  approaches  and 
the  grazing  angle  changes.  Thus,  an  integration  of  the  Illumination 
over  the  rough  surface  is  necessary.  For  example,  Beckmann  [15]  used 
isotropic  spherical  wave  Illumination  to  derive  an  expression  for  the 
coherent  field,  which  is  a confluent  hypergeometr Ic  function.  Additional 
observations  are  in  Section  6.3.  Further  work  on  this  theoretical 
problem  would  be  useful  in  experimental  work  and  for  application  to 
practical  systems. 

A fifth  factor  is  that  the  one-wav  transmitted  field 
strength  has  to  be  corrected  for  the  1/R  geometrical  behavior. 

4.3.3  Limited  Number  of  Minima 


Because  of  the  relatively  small  number  of  maxima  and 
minima  in  each  aircraft  run  (at  L-band  with  a 50-foot  receiver  height), 
it  is  difficult  to  obtain  an  adequate  statistical  sample  at  any  one  set 
of  conditions  (grazing  angle,  sea  state,  o , angle  between  the  electro- 
magnetic wave  propagation  path  and  the  ocean  wave  propagation 
direction(s) ) . Conceivably,  some  large  number  (e.g.,  ■*  100)  of  aircraft 
flights  could  be  made  in  succession  and  treat  the  "*  100  interference 
patterns  as  an  ensemble  to  obtain  the  coherent  field  and  the  incoherent 
field  using  the  statistical  signal  model  [4].  There  must  not,  however, 
be  the  large  ship-induced  oscillations  in  the  data,  described  in 
Sec.  4.2,  if  the  Incoherent  field  Is  to  be  obtained.  Even  so,  these 
> 100  flights  vrould  have  to  be  repeated  for  the  number  of  wave  heights 
desired  (o  ) and  for  a number  of  different  path  angles  relative  to 
ocean  wave'^dlr actions — an  obviously  impractical  total  number. 

In  future  experiments,  the  number  of  maxima  and  minima 
could  be  increased  as  follows.  As  the  angular  spacing  of  the  lobes  is 
approximately  \/22.,  either  the  terminal  height  (-j^)  ‘^^n  be  increased 
(e.g.,  from  the  present  50  feet  to  100  feet),  or  a shorter  wavelength 
can  be  used.  Or,  for  an  alternate  form  of  processing,  see  Sec,  6.6. 


4.3.4  Smoothed  "Apparent  Reflection  Coefficient" 


In  order  to  do  the  smoothing  of  the  interference 
patterns  to  calculate  the  coherent  field,  the  runs  are  divided  into 
equal  Intervals  of  grazing  angle.  A 0.5®  grazing  angle  Interval  con- 
tains only  one  minimum  or  maximum  at  the  lower  grazing  angles.  As  a 
compromise,  a 1"  smoothing  angle  interval  is  adopted;  a larger  Interval 
is  undesirable  because  the  coherent  field  changes  over  the  grazing 
angle  interval  with  (a  sini^)/X  , by  the  C(smooth)  vs  tp  relation,  and  as 
the  result  of  illumina¥ion  varying  with  As  C makes  up  only  part  of 
the  mean  field  at  a minimum  or  a maximum,  corrections  for  these  varia- 
tions of  C with  are  difficult.  On  the  other  hand,  the  1*  interval 
would  often  contain  only  two  minima  or  maxima.  To  obtain  a better 
statistical  sample,  given  the  experimental  setup  at  L-band  as  described, 
those  values  of  reflection  coefficient  in  the  same  grazing  angle 
interval  can  be  averaged  from  different  runs,  if  the  ocean  roughness 
and  the  propagation  path  angle  with  respect  to  the  ocean  wave  directions 
remain  constant.  This  restriction  is  not  always  satisfied  from  one  run 
to  the  next,  however,  because  either  a or  the  path  azimuth  changes. 
Internal  consistency  was  helpful  in  some  cases  as  an  aid,  as  will  be 
seen  later. 


The  signal  amplitudes  in  the  1®  Intervals  are  corrected 
for  the  1/R  factor  before  averaging.  From  the  smoothed  minima  and 
maxima  (averaged  over  the  1®  grazing  angle  intervals),  smoothed 
"apparent  reflection  coefficients"  are  calculated  from  eq.  (1).  These 
are  the  points  denoted  by  the  large  symbols  in  Figs.  5 to  9;  they 
appear  visually  to  follow  the  mean  trend  of  the  "instantaneous 
reflection  coefficient"  in  spite  of  its  large  swings. 

4.3.5  Correction  for  Presence  of  Incoherent  Field 


These  are  not,  however,  the  true  values  of  the  coherent 
field  reflection  coefficient,  because  a correction  for  the  incoherent 
field  has  yet  to  be  made.  At  an  Interference  minimum,  for  example,  the 

signal  power  measured  is  T^2  • (^D-C)*"  + T^,  where  I"  ■ incoherent  power 

and  T-  - total  power.  Because  T “■  is  larger  than  it  would  be  without 
the  incoherent  power,  the  calculated  apparent  reflection  coefficient 
will  be  smaller  than  the  true  reflection  coefficient  (C/D).  A correction 

at  a maximum,  T ‘ - (D+-C)"  + I",  is  also  required,  but  its  percentage 
change  is  small^  in  comparison  to  the  percentage  correction  at  a minimum. 

It  is  necessary  to  know,  or  assume,  the  statistical 
distribution  of  the  incoherent  field.  Rayleigh  statistics  may  be  ade- 
quate for  approximate  calculations  [4].  Where  the  sphericities  of  the 
Illuminating  and  receiving  wave  fronts  are  important,  however,  the 
incoherent  field  statistics  depart  from  a Rayleigh  distribution  [5]. 

The  procedure  is  as  follows  for  the  Rayleigh  distribution:  Calculate 

the  mean  and  the  standard  deviation  (T  ,)  of  the  minima  signal 

' m'  so  m 
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levels  in  a 1*  grazing  angle  interval.  From  their  ratio,  (T^^j)  m/Tn,. 
enter  tables  of  Rice  "signal-plus  noise"  statistics  [16]  to  obtain  the 
ratio  S/T^.  Multiply  this  ratio  by  T to  obtain  the  "signal"  phasor  at 
a minimum,  t.e.,  ” D-C.®  Perform  the  similar  calculation 

to  obtain  the  "signal" at  a®maximum,  which  is  S.  ■ (S/T  )T  ■ D+C.  The 
coherent  field  reflection  coefficient  (C/D)  is”then: 
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Although  this  procedure  is  applied  to  the  Chesapeake 
Bay  data  [1]  (f ixed-polnt-to-airplane  data),  it  can  not  be  applied  to 
the  present  ship-to-alrplane  data  (Missions  13-15)  because  of  the  non- 
random  (probably  ship  motion)  behavior  of  the  data  (Figs.  5 to  7). 

Such  calculation^  are  made  anyhow  as  a test;  in  a large  number  of  cases, 
the  ratio  (T  .)/T  of  the  minima  in  a 1*  grazing  angle  Interval  is 
greater  than  0.5227,  which  is  the  largest  value  allowable  for  Rice 
statistics  (corresponding  to  zero  coherent  field  plus  a Rayleigh 
incoherent  field).  This  result  is  consistent  with  the  observed  non- 
random  effects  in  the  data  (Figs.  5 to  7)  and  confirms  that  this 
correction  procedure  cannot  be  used  for  the  shipboard  missions. 

This  result  forces  the  use  of  the  smoothed  "apparent 
reflection  coefficients”  (Sec.  4.3.4)  as  approximations  to  the  coherent 
field  reflection  coefficients  for  the  remainder  of  this  report.  These 
are  still  adequate,  nonetheless,  for  even  if  proper  corrections  (for 
the  incoherent  field)  could  be  made,  they  would  usually  amount  to  only 
a few  percent  in  the  coherent  field  reflection  coefficient. 

5.  REFLECTION  COEFFICIENT  RESULTS 


This  section  contains  the  main  data  results  of  the  report,  the 
coherent  field  reflection  coefficient,  as  approximated  by  the  smoothed 
"apparent  reflection  coefficient"  whose  calculation  is  described  in 
Section  4.3.  Some  data  will  be  shown  first  as  a function  of  apparent 
surface  roughness  to  compare  with  Ament’s  theory  [13].  Most  of  the 
data  will  be  displayed  as  a function  of  grazing  angle  with  wave  heights 
as  parameters  to  find  how  well  wave  height  is  determined  for  remote 
sensing  purposes. 

5.1  Reflection  Coefficient  vs  Apparent  Surface  Roughness 

Figure  13  is  a comparison  of  the  reflection  coefficient  data 

of  Run  S of  Mission  15  (on  3/4/74)  with  Ament’s  theory  [13]  of  the 

coherent  field  reflection  coefficient  vs  apparent  surface  roughness. 

(o,3inv)/^.  For  the  open  square  s>'mbols,  the  values  of  (o^sin; ) .' \ are 

calculated  using  a o of  0.64  feet,  as  determined  by  the  airborne  laser 

w 
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prof llomet«r,  and  the  values  of  at  the  midpoints  of  the  1®  grazing 
angle  intervals.  The  open  squares  obviously  result  in  values  of 
abscissa  which  are  too  high;  the  solid  squares  using  j - 0.55  ft 
result  in  a good  fit  to  the  theoretical  curve.  In  this  and  all  other 
curves,  increments  of  0.05  ft  in  c are  used  as  the  smallest  warranted 
by  the  data. 

In  Run  11  of  Mission  17  (3/8/74),  a value  of  o of  0.45  ft 
provides  a good  fit  to  the  theoretical  curve,  as  shown  Yn  Fig.  14.  For 
Run  12  of  Mission  17  in  Fig.  15,  the  abscissa  values  are  calculated  with 
o ■ 0.5  ft.  In  Fig.  15,  the  aircraft  altitude  is  2500  ft;  in  all  other 
figures,  unless  speclflcallv  stated  otherwise,  the  aircraft  altitude  is 
1000  ft. 


At  the  times  of  the  runs  in  both  Figs.  14  and  15  there  are  no 

laser  profilometer  runs  to  provide  c . Estimated  values  of  o , however, 

can  be  derived  from  laser  profilometer  runs  earlier  in  the  day  as 

follows.  The  available  values  of  o„  end  about  one  hour  before  Runs  11 

w 

and  12  on  3/8/74,  as  shown  in  Fig.  16.  Wind  speed  data  exhibit  a linear 
decrease  over  the  2 1/2-hour  span  of  Mission  17  (Fig.  lb).  It  therefore 
seems  reasonable  to  extrapolate  the  c ^ data  in  a linear  trend  parallel 
to  that  of  the  wind  speed,  as  done  in  the  figure.  Some  Justification 
for  this  extrapolation  is  found  in  the  visual  sea  state  observations 
made  from  shipboard  (Appendix  C) . The  waves  are  estimated  to  have 
dropped  from  "2  to  3 ft"  at  1437  CMT  to  "1  to  2 ft"  at  1705  GMT,  or  from 
an  average  of  2.5  ft  to  an  average  of  1.5  ft.  This  40t  drop  is  consistent 
with  the  extrapolation  of  o shown  in  Fig.  lb.  The  ex'trapolated  values 
are  o 5 0.45  ft  at  the  time  of  Run  11  and  o ?:  0.40  ft  at  the  time  of 
Run  12  (to  the  nearest  0.05  ft);  these  are  to  be  compared  to  the  values 
used  for  the  abscissa  in  Figs.  14  and  15  of  0.45  ft  and  0.50  ft,  respec- 
tively. The  correspondence  is  within  the  Inaccuracies  involved  in  the 
extrapolation. 

Within  the  tolerances  shown.  Figs.  13  to  15  indicate  that  these 
particular  data  runs  are  following  Ament's  theoretical  curve  [13]  for 
roughnesses  (o  3ln.*/)/3  of  less  than  0.1.  This  is  reassuring  in  view  of 
the  possible  sKlp-caused  data  excursions  Illustrated  in  Fig.  7. 

5.2  Reflection  Coefficient  vs  Grazing  Angle  Plots  for  Remote  Sensing 
of  Wave  Height 

IJhen  the  measured  reflection  coefficients  are  plotted  on  graphs 
of  Ament's  [13]  theoretical  curves  as  a function  of  grazing  angle  with 
wave  height  (o  ) as  a parameter  (Fig.  17),  the  values  of  to  secure  a 
fit  are  automatically  displayed.  This  results  in  an  easy  graphical 
method  of  finding  wave  heights  for  remote  sensing.  A series  of  these 
curves  (with  \ corresponding  to  a frequency  of  1.3  GHz)  follows  to  test 
how  well  the  wave  heights  determined  by  this  method  agree  with  wave 
heights  given  by  the  laser  prof ilometer . Table  Xi’  will  summarize  the 
results  of  this  series  of  figures. 
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The  data  of  Figs.  13  and  13  are  replocted  in  this  format  for 
comparison.  Figure  18  (Run  8,  Mission  15)  is  one  of  the  best  data  runs 
in  terms  of  length  of  run  and  the  small  scatter  of  points  using  optical 
tracking.  The  predicted  wave  height  (o  ) is  easily  determined  by  the 
points  as  o ■ 0.55  + 0.05  ft,  compared'^to  0.64  ft  from  the  laser 
prof ilometer . (Again,  increments  of  0.05  ft  are  the  smallest  which 
will  be  used.) 

Similarly,  Fig.  19  (Run  11,  Mission  17)  is  another  example  of 
a good  run  in  which  a can  be  determined  by  cursory  inspection  as 
o ■ 0.45  + 0,05  ft.  (The  point  at  a grazing  angle  of  1*  is  obviously 
in  error.)  The  extrapolated  profilometer  a is  0.45  ft  (Sec.  5.1). 

V 

In  Fig.  20  (Run  12,  Mission  17)  with  only  four  points,  the 
scatter  is  wider,  but  o is  determined  as  J * 0.5  ft  within  the  toler- 
ance of  0.05.  The  extrapolated  profilometer  c is  0.40  ft  (Sec.  5.1). 
— ^ 

Thus,  in  these  three  runs,  the  reflection  coefficient  data 
determined  the  wave  height  standard  deviation  (o  ) to  a tolerance  of 
+ 0.1  ft. 


Figure  21  (Run  5,  Mission  15)  shows  considerable  scatter,  but 
two  points  closely  follow  the  o ■ 0.6-ft  curve  and  two  points  fail 
around  the  o • 1.0-  to  1.1-ft  curves.  A possible  cause  of  this 
behavior  migKt  be  how  the  large  data  excursions  (such  as  those  which 
appear  in  the  data  obtained  by  optical  tracking  in  Mission  15)  happen 
to  straddle  the  1*  smoothing  intervals.  In  Run  8 of  Mission  15  (Figs.  7 
and  18),  apparently  the  excursions  meshed  with  the  1*  intervals  in  such 
a way  as  to  almost  average  them  out.  In  general,  a bias  upward  in  some 
1®  intervals  and  downward  in  others  might  be  expected,  given  the  large, 
regular  excursions  of  the  optical  tracking  data  (Fig.  7).  The  average 
0 of  the  five  points  in  Fig.  21  is  0,8  + 0.3  ft,  compared  to  an 
interpolated  laser  profilometer  o of  0.8  ft.  The  exact  agreement  is 
fortuitious,  of  course. 

The  next  three  Figures  (22,  23,  and  24)  are  data  from  Mission  16 

on  3/6/74.  Although  there  are  no  laser  profilometer  data  on  this  day, 

visual  sea  state  observations  made  from  shipboard  (.\ppendix  C)  are 

"swells  approximately  3 ft,  occasionally  5 ft."  Assuming  that  the  3-ft 

value  corresponds  to  the  significant  wave  height,  then 

••  ft-  Similarly,  assuming  that  “ 5 ft,  then 

0^  ■ (Hj^yj^p) /5. 1 1 ft.  Thus,  visual  observations  give  a range  for 

of  approximagely  O.T5  to  1 ft.  In  Fig.  22,  the  points  of  Run  2 

indicate  an  average  value  of  0 , of  0.65  + 0.15  ft,  or  0.1  ft  below 

the  side  of  the  visual  ranges  of  values.  In  Fig.  23,  Run  5 yields  a 

0 of  0.9  + 0.1  ft,  which  falls  within  the  visual  range  of  values, 
w — 

In  Fig.  24  are  the  remaining  two  Mission  16  data  runs  which  do 
not  have  any  internal  consistency  and  are  widely  scattered.  The 
values  of  indicated  for  the  Run  3 and  4 points  are  1.0  and  0.6  ft. 


respectively  (about  the  same  spread  given  by  the  other  Mission  In  runs 
in  Figs.  12  and  23),  and  the  average  for  all  the  points  in  Fig.  24  is 
■ 0.3  + 0,2  ft.  Thus,  the  average  o ^ derived  from  all  the  Mission 
16  points  is  0.8  + 0.2  ft,  which  is  at  Yeast  consistent  with  the 
visually  estimated  values  of  0.75  to  1.0  ft.  If  the  actual  value 
happened  to  be  0.8  ft,  the  error  of  + 0.2  ft  would  correspond  to  + 25t. 

The  tracking  radar  was  repaired  and  placed  into  operation  for 
the  first  time  during  Mission  16.  Run  2 used  optical  tracking  to  point 
the  receiving  horns,  but  radar  tracking  pointed  the  horns  for  Run  3 and 
thereafter.  It  is  Interesting  to  note  the  very  wide  scatter  of  points 
in  Run  3,  medium  scatter  in  Run  4,  and  the  low  scatter  in  Run  5.  Is 
this  a result  of  improvement  in  radar  performance,  after  Just  being 
turned  on  at  1243  GMT,  from  the  beginning  of  Run  3 at  1252  GMT  to  the 
end  of  Run  5 at  1317  GMT?  Or  just  happenstance?  .At  least  in  Mission 
18,  with  a different  radar  on  a different  ship,  these  problems  of  wide 
scatter  are  relatively  minor.  This  mission  will  be  discussed  next. 

In  contrast  to  Fig.  24,  the  much  reduced  scatter  of  points 
with  good  radar  tracking  is  shown  in  Fig.  25  (Runs  1 and  3,  Mission  18). 
In  Run  1,  one  of  the  four  points  departs  considerably  from  the  trend  of 
the  other  three;  it  is  very  similar  in  appearance  to  the  radar  tracking 
run  in  Fig.  20.  In  Run  3 (Fig.  25),  only  one  point  (at  * ■ 1.'>*1  of 
the  five  is  obviously  in  error;  this  run  is  similar  to  the  radar 
tracking  run  in  Fig.  19,  which  also  contains  a low  point  at  the  same 
angle.  Thus,  the  radar  tracking  runs  may  have  perhaps  one  point 
clearly  in  error,  but  the  remainder  of  the  points  form  a consistent 
trend  with  relatively  small  scatter.  The  points  of  Run  1 give  an 
average  ■ 0.06  + 0.15  ft;  the  points  of  Run  3 ^excluding  the  point 
at  »’  " 1.5*)  yield  ■ 0.05  + 0.1  ft.  The  average  value  of  the  four 

laser  profilometer  values  (0.65,  0.72,  0.5*^,  and  0.o2  ft)  taken 

during  the  period  of  Runs  1 to  4 is  ■ 0.o5  + 0.07  ft,  which  agrees 
within  the  error  tolerances. 

The  roughest  twv  days  are  for  Missions  14  and  13.  There  is 
not  much  to  say  about  Mission  14  in  Fig.  2o  except  that  it  is  a bad 
example.  Although  there  are  six  runs,  the  number  of  points  per  run 

ranges  from  only  1 to  3,  not  enough  to  use  internal  consistency  as  an 

aid.  Even  the  averages  of  the  points  at  each  grazing  angle  do  not  have 
the  correct  trend.  This  is  an  optical  tracking  day,  and  the  large 
scatter  is  presumed  to  be  a result  of  the  large  pointing  excursions,  as 
previously  explained.  The  average  o given  by  all  of  the  points  in 
Fig.  2n  is  ^ 0.3  + 0.5  ft,  compared  ¥o  the  average  laser  profilometer 
value  of  1.15  + 0.12  ft.  or  30*  low.  a large  error. 

The  data  for  the  roughest  day,  2,'25/''-*  (Mission  13),  are  given 
in  Fig.  27.  The  wide  scatter  follows  from  the  optical-tracking  behavior 
of  the  instantaneous  values  shown  in  Figs.  5 and  n.  There  are  only  two 
points  in  two  runs  and  three  points  in  the  third  run;  a run  is  not  long 
enough  to  use  Internal  consistency  as  a help.  The  indicated  values  of 
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range  from  0.7  to  1.2  ft,  but  with  no  proportional  correspondence  to 
the  profilometer  values,  i.e,,  a behavior  more  consistent  with  Just 
random  errors.  The  average  of  the  values  from  all  the  points  in 
Fig.  27  is  0,9  + 0.3  ft,  compared  to  an  average  of  1.33  ft  for  the 
profilometer  values,  or  a large  error  of  32*. 

The  results  of  the  comparisons  made  in  Figs.  18  to  27  are 
listed  in  detail  in  Table  K\’,  including  the  dubious  runs.  This 
comparison  will  be  discussed  next  in  Section  6. 

6.  DISCUSSION 

6.1  Errors  in  the  Remotelv  Sensed  Values  of  o.. 

: — — H. 

In  assessing  the  validity  of  using  reflection  coefficients  as 
a remote  sensing  indicator  of  cr^,,  it  is  necessary  to  use  data  runs  in 
which  there  is  most  confidence.  These  "qualified"  runs  have  4 or  more 
points  w'ith  some  internal  consistency;  they  are  summarized  in  the 
shorter  Table  XVI  for  ease  in  comparison.  The  average  percent  error 
for  all  runs  with  4 or  more  points  per  run  just  happens  to  balance  out 
to  zero,  but  could  have  been  + 5*.*  In  Mission  16  when  only  visual 
estimates  of  wave  heights  were  available,  an  error  estimate  of  less  than 
+ 25*  was  made;  this  value  was  not  Included  in  the  averages  because 
neither  the  magnitude  nor  the  sign  of  the  error  was  known.  If,  however, 
the  Mission  lb  error  were  the  maximum  of  +255,  the  average  error  for 
the  runs  vrould  increase  from  05  to  +35;  if  the  Mission  16  error  were 
-255,  the  average  error  for  the  runs  would  change  from  05  to  -45, 
neither  change  being  very  significant. 

The  error  limits  on  the  deduced  values  of  c^,  (third 
column.  Table  X\’I3  encompass  the  maximum  spread  of  the  points.  For  the 
five  radar-tracked  runs,  these  are  + 0.15  ft  for  one  run  and  + ^ 0.1  ft 
for  the  remaining  four  runs.  Even  these  moderate  tolerances  result  in 
a 255  error  at  the  lowest  wave  height  of  0.4  ft  (.Run  12  of  Mission  17'), 
whereas,  at  the  largest  of  O.Q  ft  for  radar-tracked  runs  (Run  5, 
Mission  16),  the  ^ 0.1-ft  tolerance  causes  only  an  115  error. 

Thus,  the  comparison  of  obtained  from  the  "qualified" 
reflection  coefficient  runs  to  the  laser  profilometer  c^.'s  is  favorable. 
Even  if  the  optically-tracked  runs  of  Missions  13  and  14,  which  are 
short  and  of  doubtful  validity,  are  thrown  in  for  a worst  cast  estimate, 
the  average  error  for  all  runs  is  still  -145  (right-hand  column  of 
Table  KVl.** 


* Taking  the  worst  cases  in  Column  5 of  Table  X\'I  and  assuming  that 
either  only  the  +255  error  existed  for  the  6 runs,  or  that  only  the 
-225  error  existed. 

**  This  is  the  result  of  assuming  that  the  total  of  21  points  in  Missions 
13  and  14  is  the  equivalent  of  5 runs  of  4 points  each.  These  five 
runs  with  an  error  of  -315  per  run  are  then  weighted  with  the  remainder 
of  the  runs  in  Table  X\’I  to  obtain  -145. 
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On  the  two  roughest  days,  Missions  13  and  14,  the  reflection 
coefficient  data  indicate  values  of  which  are  approximately  31%  low, 
i.e.,  the  reflection  coefficients  are  higher  than  they  should  be  for 
correspondence.  Is  this  the  result  of  a trend  at  higher  surface 
roughnesses,  or  a result  of  the  large,  regular  data  excursions  resulting 
(presumably)  from  optical  tracking?  The  highest  profilometer  value  of 
on  any  other  day  (for  which  reflection  coefficients  exist)  is  0.8  ft 
from  Run  5 of  Mission  15  (an  optical  tracking  day);  the  reflection 
coefficient  data  yield  • 0.8  ft,  to  agree  with  the  profilometer 
value.  On  Mission  16,  another  day  of  comparable  roughness,  visual  wave 
observations  give  as  0.75  to  1 ft;  the  average  of  all  values  from 
the  Mission  16  reflection  coefficients  is  0.8  + 0.2  ft,  a consistent 
value.  Thus,  any  trend,  if  one  exists,  would  have  to  begin  for  values 
of  Oy.  above  0.8  ft  and  would  have  to  be  30%  low  by  a of  1.14  ft 
(Mission  14).  Such  a trend  seems  doubtful;  the  values  of  roughness 
(a^fSini;/)/!  reached  in  the  data  in  this  report  are  generally  4 0.1,  and 
the  earlier  data  taken  between  fixed  points  [4]  (where  ship  and  airplane 
motions  are  not  involved)  show  an  experimental  trend  of  reflection 
coefficient  breaking  away  above  Ament's  [13]  theoretical  curve  only  for 
(a^sinil') /A>0. 1.  Thus  it  seems  more  likely  that  the  30%  low  values  of 
Oy  given  by  Missions  13  and  14  are  the  result  of  the  large,  optical 
tracking-induced  excursions.  Moreover,  it  would  be  expected  that  the 
most  difficulty  with  optical  tracking  would  occur  on  the  roughest  days. 

6.  2 Methods  of  Reducing  Error  in 

There  are  several  feasible  methods  of  reducing  the  error  in 
the  individual  runs  of  the  method  so  as  to  approach  the  low  average 
error  that  is  realized  for  all  the  runs.  The  major  error  source  is  an 
insufficient  number  of  minima  and  maxima  within  the  grazing  angle 
smoothing  interval. 

a.  Operate  at  a shorter  wavelength,  such  as  3 cm.  This  would 
produce  seven  times  as  many  minima  in  the  same  grazing  angle  interval 
as  the  1.3  GHz  used  in  these  tests.  (A  new  set  of  Ament's  curves  has 
to  be  drawn  for  the  new  wavelength.) 

b.  Increase  the  size  of  the  grazing  angle  smoothing  interval 
to,  say,  2"  instead  of  the  1"  used  in  this  report.  This  change  makes 
the  corrections  listed  in  Section  4.3.2  more  difficult.  A new  set  of 
Ament's  curves  vs  grazing  angle,  but  averaged  over  2®  intervals,  also 
is  required. 


c.  Raising  the  lower  terminal  height  (the  shipboard  antenna 
herein)  will  also  increase  the  number  of  minima  in  the  grazing  angle 
interval. 


d.  Use  wide  antenna  beamwidths  in  any  case  to  avoid  problems 
of  pointing  errors. 
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e.  For  best  accuracy,  choose  a range  of  grazing  angles  such 
that  there  is  a wide  spacing  between  the  different  a curves  over  the 
range  of  wave  heights  to  be  covered  (such  as  4®  to  5®  in  Fig.  17  for 
1.3  GHz),  but  such  that  the  minimum  reflection  coefficient  is  greater 
than  0.4.  The  reason  for  the  latter  stipulation  is  that  Ament's  [13] 
theoretical  curves  hold  for  Ccr„sint<))/A^0.1,  but  the  experimental  data 
"curve"  [4]  breaks  away  and  becomes  larger  than  Ament's  theoretical 
curve  for  Ca^slnii;)/X>0.1.  This  value  of  roughness  corresponds  to  a 
coherent  field  reflection  coefficient  of  0.4.  As  there  are  conflicting 
theories  for  the  coherent  field  behavior  above  roughnesses  of  0.1,  it  is 
necessary  at  the  moment  for  remote  sensing  applications  to  restrict  the 
measured  reflection  coefficients  to  ^ 0.4  so  that  Ament's  [13]  theoreti- 
cal curves  can  be  used. 

6.3  Non-Uniform  Illumination  Errors 

The  receiving  horn  elevation  beamwidths  of  32®  were  sufficiently 
wide  when  the  beams  were  pointed  horizontally  as  in  the  Chesapeake  Bay 
tests  [1].  On  shipboard,  the  horns  were  mounted  on  top  of  the  MK  25 
tracking  radar  antenna  to  obtain  a mounting  point  that  was  cross-level 
stabilized.  In  tracking  the  aircraft,  the  radar  pointed  the  horns  at 
the  aircraft,  but  also  elevated  the  beams  and  reduced  the  illumination 
of  the  water  at  the  higher  grazing  angles.  There  is  no  simple  way  to 
correct  for  this  effect  (Sec.  4.3.2).  Some  authors  have  made  a multi- 
plying "correction"  factor  from  relative  strengths  of  the  direct  ray  and 
the  specularly  reflected  ray  directions  as  they  emerge  from  the  antenna 
patterns,  but  this  procedure  is  not  necessarily  correct  [4]. 

In  Fig.  28,  the  points  denoted  by  open  squares  are  the  same  as 
given  before  in  Fig.  18.  These  points,  after  being  multiplied  by  the 
ratio  of  direct  to  specular  ray  antenna  beam  amplitudes,  for  both  trans- 
mitting and  receiving  beams,  become  the  asterisks.  The  two  asterisks  at 
the  two  highest  grazing  angles  appear  to  be  "overcorrected,"  and  the 
average  o^  indicated  by  the  asterisks  is  0.45  ft,  compared  to  0.55  ft 
for  the  uncorrected  points  and  0.64  ft  for  the  laser  prof ilometer . In 
Fig.  29,  the  "overcorrections"  are  obvious,  with  the  asterisks  indicating 
a 0^  of  0.25  ft,  compared  to  the  of  0.45  ft  given  by  both  the  uncor- 
rected points  and  the  laser  prof ilometer . Thus,  the  uncorrected  points 
agree  with  the  prof ilometer  far  better  than  the  so-called  "corrected" 
points.  There  is  a possible  reason  for  part  of  this  behavior.  Since 
these  receiving  horns  were  aligned  visually,  there  is  at  least  the 
possibility  that  they  were  pointed  a few  degrees  downward  instead  of  at 
zero,  as  reported,  when  the  radar  dish  was  pointed  at  zero  degrees 
elevation.  This  would  mean  that  the  illumination  was  greater  than 
thought,  and  that  any  "corrections"  would  be  too  much,  as  Figs.  28  and 
29  Indicate.  Thus,  the  reflection  coefficients  in  all  the  figures  prior 
to  Fig.  28  are  uncorrected  ones,  as  measured. 
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This  problem  o£  the  quantitative  effect  of  non-uniform  Illumi- 
nation on  the  rough  surface  coherent  field  Is  one  requiring  theoretical 
work.  Meanwhile,  the  experimenter  needs  to  use  wide  beams  to  avoid  the 
effect. 


6.4  Directional  Effect 


There  are  Insufficient  data  to  find  wliether  there  is  any 
directional  effect,  l.e.,  of  the  microwave  propagation  path  relative 
to  the  ocean  wave  propagation  direction(s) . Of  the  runs  of  4 or  more 
points  In  Table  X\'I,  In  only  two  is  the  aircraft  heading  relative  to 
the  wind  known:  Run  8,  Mission  15,  Is  upwind  at  75*  to  the  wind 
(in\r/75*),  and  Run  11,  Mission  17,  is  downwind  at  30“  to  the  wind 
(DW/30*).  In  both  of  these  runs,  the  predicted  values  of  agreed 
with  the  laser  profllometer  values  within  +0.1  ft.  The  otlier  short 
runs  of  Missions  13,  14,  and  15  (listed  in  Table  XV)  give  contrad Ictorv 
results,  l.e.,  there  are  more  errors  in  these  short,  optically-tracked 
runs  than  whatever  directional  effect  might  exist.  Thus,  a large 
number  of  flights  would  be  necessary  to  establlsli  and  measure  a direc- 
tional effect  by  this  method. 

6.5  Laser  Profllometer  Errors 


The  Information  on  the  laser  profllometer  performance  is  as 
follows.  Two  Missions  (14  and  18)  provided  a good  test  of  repeata- 
bility In  that  in  each  of  the  Missions,  the  values  of  from  the 
different  runs  were  tightly  clustered  and  sliowed  no  trend  (Appendix  B). 
The  entries  of  profllometer  In  Tables  X\’  and  XVT  for  Missions  14  and 
18  are  the  means,  plus  or  minus  the  spread  of  all  the  values  about  the 
mean  for  all  of  the  runs  (on  each  day).  If  the  profllometer  is  assumed 
(for  discussion)  to  have  no  error,  the  small  spread  shows  the  constancy 
of  the  waves;  conversely,  if  the  sea  is  assumed  constant,  the  spread 
gives  the  repeatability  of  the  laser  profllometer  measurements.  For 
both  days  the  repeatability  spread  is  + 10*. 

Information  on  the  absolute  accuracy  of  the  profllometer  is 
scanty.  Measurements  by  a laser  profllometer  have  been  compared  with 
those  of  a wave  staff  by  other  organlrat ions  with  satisfactory  results. 
Although  photographs  of  the  ocean  taken  from  the  airborne  L-band 
transmitter  pod  were  submitted  to  the  appropriate  group  for  spectral 
analysis,  the  report  received  was  that  no  wave  spectra  could  be  obtained 
for  a variety  of  reasons,  such  as  sun  glitter,  too  great  an  angle 
between  wind  and  wave  directions,  not  enough  waves  In  the  photograph, 
etc.  In  these  tests,  the  course  of  the  S-2D  aircraft  carrying  the 
transmitter-camera  pod  was  constrained  by  tlie  slilp's  course,  and, 
consequently,  photographs  may  not  have  been  taken  under  ideal  condi- 
tions. At  least  In  these  tests,  the  requirements  for  taking  photographs 
for  analysis  bv  the  Stllwell  method  were  too  constraining  and  restric- 
tive to  allow  results  to  be  obtained. 


A comparison  of  profilometer  results  to  visual  observations 
(Appendix  C)  shows  the  profilometer  values  being  about  half  as  large 
iMission  13)  to  twice  as  large  (Mission  18)  with  reasonable  agreement 
in  between  (Mission  17),  This  is  probably  consistent  with  the  accuracy 
of  the  visual  observations. 


A circuitous  reasoning  provides  some  confidence  in  the  laser 

profilometer  data,  however.  As  the  reflection  coefficients  agree  with 

Ament's  curve  when  the  profilometer  o^.'s  are  used  to  calculate  apparent 

surface  roughness,  and  because  Ament's  curve  has  been  confirmed  in  the 

past,  then  the  profilometer  values  of  o should  be  correct. 

w 

Laser  profilometer  spectra  could  not  be  used,  however,  because 
of  aliasing  inherent  in  making  a line  measurement  of  a two-dimensional 
surface  spectrum.*  This  does  not  affect  the  values  of  o . 

6,6  Proposed  FM-Sweep  Method  and  Remote  Sensing  Application 

In  the  data  processing  used  herein,  a statistically  significant 
number  of  maxima  and  minima  can  not  be  obtained  In  the  l'’-gra2ing-angle 
intervals  without  an  inordinate  number  of  aircraft  flights  (Secs.  4.3.3 
and  4.3.4). 


To  avoid  these  difficulties,  an  FM-sweep  technique  is  proposed, 
in  which  the  microwave  frequency  is  swept  at  a rapid  rate,  thus 
presenting  (on  an  oscilloscope,  e.g.)  an  essentlallv  "instantaneous" 
section  of  the  interference  pattern.  The  FM-sweep  rate  is  set  high 
enough  so  that  signal  changes  do  not  occur  within  one  sweep  period,  and 
thus  the  waves  can  be  considered  to  be  "frozen."  Then  a sufficient 
number  of  such  sweeps,  depending  upon  geometrv  and  velocitv.  is  recorded 
to  obtain  an  adequate  ensemble  for  statistical  treatment. 


The  signal  processing  is  straightforward.  The  frequency 
deviation  of  the  FM  can  be  adjusted  to  sweep  through  one  maximum  and  one 
minimum  of  the  Interference  pattern  on  eacli  sweep.  One  could  imagine 
many  sweeps  being  integrated  by  a long-persistence  oscilloscope  screen, 
with  photographs  being  taken  of  the  averaged  interference  pattern  traces, 
but  this  is  an  unwieldy  method.  A possible  analog  processing  method 
would  be  to  use  two  envelope  detectors,  back-to-back,  the  bias  of  one  set 
for  detection  of  the  signal  maxima,  and  the  bias  of  the  reverse  detector 
set  for  detection  of  the  signal  minima.  For  remote  sensing,  a simple 
procedure  is  to  pa^s  the  two.  detector  outputs  through  long-time  constant 
filters  to  obtain  ^mln*  T’hen,  C/P  is  calculated  from  the 

approximation: 


C/D  * (T, 


Max 


^min/' 


iT, 


Max 


^min' 


* Private  communication  from  Donald  L.  Hammond,  NRL. 


Plotting  C/D  on  a set  of  curves  like  Fig.  17  (but  calculated  for  the 
value  of  \ being  used)  determines  of  the  waves.  A given  value  of 
grazing  angle  (and  A)  determines  a vertical  slice  through  the  curves 
In  Fig.  17.  An  equation  can  be  determined  for  this  slice,  so  that  the 
value  of  C/D  could  yield  directly  without  plotting,  if  desired. 

The  antenna  beamwldth  needs  to  be  constant  over  the  frequency 
sweep  (e.g.,  10  to  12  GHz).  The  system  gain  should  preferably  also  be 
constant  over  this  frequency  band,  although  a calibration  curve  correc- 
tion could  be  inserted  in  the  computations. 

This  FM-sweep  method  trades  some  complexity  of  equipment  and 
data  handling  for  a prohibitive  number  of  aircraft  flights.  Even  though 
this  trade  would  usually  be  desirable,  the  FM  method  in  addition  has  the 
unique  feature  of  presenting  "instantaneous"  sections  of  Interference 
patterns  for  "frozen"  seas. 

Would  it  be  suitable  for  remote  sensing  by  forward  scatter 
from  satellites,  however?  There  are  two  obvious  drawbacks.  One  is 
poor  resolution  because  the  Fresnel  zones  cover  a large  area  at  grazing 
angles  of  4“  to  6*.  Because  It  has  been  shown  that  20  Fresnel  zones 
provide  all  but  a few  percent  of  the  scattered  energy,  the  "resolution 
area"  could  be  considered  as  the  area  inside  the  20th  Fresnel  zone. 

The  second  point  is  that  although  forward  scattering  requires 
two  satellites,  one  each  for  transmitter  and  receiver,  this  is  not  the 
devastating  requirement  that  it  once  was.  A system  which  meets  the 
geometry  needed  for  the  remote  sensing  methods  in  this  report  is  being 
proposed  by  the  Applied  Physics  Laboratory  of  Johns  Hopkins  University 
[21]  to  measure  surface  pressure  for  weather  forecasting.  It  consists 
of  two  satellites  which  make  use  of  both  the  direct  (line-of-slght  path) 
between  them  and  also  the  specularly  reflected  path  from  the  earth's 
surface  at  a grazing  angle  of  4®  to  6®.  It  Is  conceivable  that  such 
satellites  might  be  designed  to  also  carry  the  equipment  for  the  sea 
surface  remote  sensing  as  proposed  in  this  report  and  thus  serve  a dual 
application. 

On  the  positive  side,  the  accuracy  in  wave  height  obtained  with 
the  method  used  herein  (l®-slices  In  grazing  angle)  compares  very  favor- 
ably with  the  preflight  specifications  of  the  altimeter  in  SEASAT-A,  and 
the  proposed  FM-sweep  method  should  increase  this  accuracy.  In  addition, 
procedures  for  possible  remote  sensing  of  the  two-dimensional  corre- 
lation lengths  of  the  ocean  surface  will  be  described  In  another  report. 


7.  SUMMARY  AND  CONCLUSIONS 


7.1  Remote  Sensing  of  Significant  Wave  Height  (4o  ) 

w 


1.  The  measured  reflection  coefficients  from  the  six  qualified* 
runs  follow  Ament’s  [13]  theoretical  curve  of  reflection  coefficient  as 

a function  of  apparent  surface  roughness  (o^^sini|;)/X.  Only  small  rough- 
nesses of  4 0.1  are  encountered  at  this  microwave  frequency  of  1.3  GHz. 

2.  As  a remote  sensing  method,  a series  of  Ament's  [13] 

theoretical  curves  is  recast  as  a function  of  grazing  angle,  with  wave 

height  (o^)  as  a parameter.  A measured  reflection  coefficient  value 

plotted  at  a given  grazing  angle  then  gives  a directly. 

w 

3.  The  average  spread  in  the  values  of  predicted  by  this 
remote  sensing  method  is  + 0.12  ft  in  the  six  runs.  This  is  equivalent 
to  errors  in  individual  runs  of  from  + 30%  at  the  lowest  wave  height 
(0.4  ft)  to  + 13%  at  the  highest  wave  height  (0.9  ft)  encountered. 

4.  The  accuracy  of  this  remote  sensing  method  is  indicated  by 
the  fact  that  the  average  (over  the  six  runs)  percentage  error  of  the 
predicted  compared  to  the  laser  profilometer  is  only  + 5%. 

5.  Several  feasible  methods  of  reducing  the  error  in  the 
individual  runs  of  the  method  so  as  to  approach  a low  overall  average 
error  are  listed.  The  major  error  source  is  an  insufficient  number  of 
minima  and  maxima  within  the  l®-grazing  angle-smoothing  interval.  The 
best  remedy  is  to  operate  at  a shorter  wavelength,  such  as  3 cm.  This 
would  produce  seven  times  as  many  minima  in  the  same  grazing  angle 
interval  as  at  the  1.3-GHz  frequency  used  in  these  tests. 

6.  There  are  insufficient  data  to  determine  whether  any 
directional  effect  exists,  i.e.,  angle  of  electromagnetic  wave  propa- 
gation path  relative  to  the  ocean  wave  direction(s) . 

7.  There  is  not  an  adequate  statistical  data  base  for  calcu- 
lating the  Incoherent  field  because  of  the  large,  regular  excursions 
in  the  optically-tracked  data,  and  because  of  an  inadequate  number  of 
minima  in  the  grazing  angle-smoothing  interval  at  L-band. 

7.2  Proposed  FM  Sweep  Method  for  Remote  Sensing  of 


A system  is  proposed  for  remote  sensing  of  o^  which  repeti- 
tively sweeps  (FM)  the  microwave  carrier  frequency  at  a rapid  rate  to 


* The  "qualified"  runs  are  six  runs  with  four  or  more  points  and  with 
internal  consistency.  These  are  primarily  only  the  radar-tracked 
runs,  but  include  two  optically-tracked  runs  which  have  internal 
consistency. 
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present  an  "instantaneous"  interference  pattern  that  "freezes"  the  waves 
as  far  as  signal  changes  are  concerned.  An  ensemble  average  is  then 
performed  to  calculate  the  reflection  coefficients  and  o„  with  only  a 
relatively  few  aircraft  flights  for  different  azimuths. 
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0.79 
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0.90 

IS 
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0.61 

17 

5.18 

0.69 

20 

4.40 

0.49 

19 

4.63 

0.53 

2'» 

4.00 

0.40 

21 

4.19 

0.44  1 

24 

3.67 

0.33 

23 

3.33 

0.36 

26 

3 39 

0.23 

25 

3.52 

0.31 

23 

3.14 

0.25 

n* 

3.26 

0.26 

30 

2.93 

0.21 

29 

3.04 

0 22 

32 

2.75 

0.18 

31 

2.34 

0.20 

34 

2.59 

0.16 

33 

2.67 

0.17 

Table  VI 

Comparison  of  Measured  and  Calculated  Minima  Positions 


Table  VII 

Phase  Changes  Produced  by 
Aircraft-Height  Variations 


AZo  s 


“ 2 Zj  N 


R 


Z^  * 51ft.  X - 0.757  ft.  (1.3  GHz> 


Range 

Altitude  Variations  AZ» 

(ft.) 

R 

to  cause  phase  changes  of 

ft. 

180'* 

90° 

45° 

O 

O 

CO 

10,000 

37 

19 

9 

6 

75,000 

278 

139 

70 

46 

Table  VIII 

Signal-Level  Changes  vs.  Phase  Angle  from  a Minimum 


5 

a 

P 

^T'^D 

Signal  Level 

Signal  Level 

Phase  angle 

Phase  angle 

Reflection 

relative  to 

relative  to  a 

from  a 

between 

Coefficient 

direct  signal 

maximum 

minimum 

direct  and 

(assumed) 

reflected 

signals 

(deg) 

(deg) 

dB 

dB 

0 

180 

0.9 

0.100 

-20 

-25.6 

15 

165 

0.26 

-11.7 

-17.3 

30 

150 

0.500 

-6.0 

-11.6 

45 

135 

0.735 

-2.7 

-8.3 

66 

114 

1.000 

0.0 

-5.6 

1.345 

+2.6 

-3.0 

180 

0 

1.900 

+5.6 

0.0 

Calculations  assume  a smooth,  plane  surface. 

Ej'Ed  * [1  P"  + 2 p cos  a]' 
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Table  X 

Interference  Pattern  Data,  Mission  13.  Run  9 


Time  - GMT 

Range 

ft 

w 

deg 

Max. 

dB 

.Min. 

dB 

Diff. 

dB 

"Instantane- 
ous Reflection 
Coefficient” 

min 

sec 

1 

13 

23 

20.10 

37807 

1.60 

— 

-82.0 

— 

29.40 

35609 

1.69 

-43.3 

— 

38.7 

0.976 

34.40 

34426 

1.75 

— 

-67.6 

24.3 

0.884 

39.10 

33315 

1.81 

-43.0 

— 

24.6 

0.888 

43.40 

32298 

1.86 

— 

-64.4 

21.4 

0.842 

44.50 

32038 

1.88 

-40.8 

— 

23.6 

0.376 

49.20 

30927 

1.95 

— 

-61.7 

20.9 

0.832 

30.70 

30572 

1.97 

-47.0 

— 

14.7 

0.690 

51.60 

30359 

1.98 

— 

-72.3 

25.8 

0.902 

33.10 

30005 

2.00 

-41.0 

— 

31.8 

0.949 

i 

56.70 

29154 

2.06 

— 

-66.7 

25.7 

0.901 

1 

59.40 

28515 

2.11 

-41.2 

— 

25.5 

0.900 

i 13 

24 

02.10 

27877 

2.16 

-60.0 

18.8 

0.792 

) 

02.80 

27711 

2.18 

-43.8 

— 

16.2 

0.732 

! 

04.40 

27333 

2.21 

— 

-62.5 

18.7 

0.790 

1 

05.60 

27049 

2.23 

-45.0 

— 

17.5 

0.775 

1 

07.60 

26576 

2.27 

— 

—50.4 

3.4 

07,80 

26529 

2.27 

—45.7 

— 

4.7 

09.90 

26033 

2.32 

— 

-63.3 

17.6 

0.776  1 

10.60 

25867 

2.33 

-44.7 

— 

19.1 

0.800 

12.20 

25489 

2.36 

— 

—67.5 

23.3 

0.872  ! 

15.20 

24779 

2.43 

—40.5 

27.0 

0.913 

1 

15.70 

24661 

2.44 

— 

-73.3 

32.8 

0.954  , 

1 

17.20 

24306 

2.48 

-47.0 

— 

26.3 

0.907  i 

1 

18.30 

23928 

2.52 

— 

-63.0 

16.0 

0.728 

1 

! 

20.10 

23621 

2.55 

—45.2 

— 

17.8 

0.770 

1 

20.60 

23503 

2.56 

-62.5 

17.3 

0.760 

I 

1 

22.90 

22959 

2.62 

-42.5 

20.0 

0.317 

' 

23.10 

22911 

2.63 

— 

—55.8 

13.3 

0,650  , 

23.95 

22710 

2.65 

-42.9* 

— 

12.- 

0.60  ; 

1 

24.80 

22510 

2.68 

— 

-61.4 

18.5 

0.786 

1 

23.20 

22415 

2.69 

-43.3 

— 

IS.l 

0,777  ‘ 

; 

26.30 

22155 

2.71 

— 

—65.3 

22.0 

0.852  j 

i 

27.10 

21966 

2.73 

-43.6 

— 

21.7 

0.847  ' 

1 

23.00 

21753 

2.77 

— 

— 67.5 

23.9 

0.880  1 

i 

29.00 

21516 

2.80 

-42.3 

— 

24.7 

0.890 

29.40 

21422 

2.81 

— 

—66.6* 

24.- 

0.88  ! 

29.80 

21327 

2.83 

-42.3 

— 

24.- 

0.88  ! 

1 

30.60 

21138 

2.85 

— 

-65.8 

23.5 

0.875 

1 

31.30 

20973 

2.86 

-40,8 

— 

25.0 

0.893 

1 

32.40 

20713 

2.90 

— 

-64.2 

23.4 

0.873 

* ■ Missina  values  "faked”  as  the  average  of  the  two  surroundina  values  to  make  the 
computer  proaram  run. 
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Table  XI 

Interference  Pattern  Data.  Mission  13.  Run  12 


Time  - GMT 

Raniie 

ft 

deg 

Max. 

dB 

— 

Min. 

dB 

Diff. 

dB 

“Instantane-  | 
ous  Reflection 
Coefficient"  | 

mm 

mm 

sec 

18 

46 

19.70 

1 

59315 

1.01 

— 

•> 

1 

58.50 

49205 

1.22 

-43.3 

— 

— 

— 

47 

15.40 

44801 

1.34 

— 

-72.5 

29.2 

0.931 

33.60 

40069 

1.50 

-45.8 

— 

26.7 

0.910 

44.00 

37349 

1.61 

-78.3 

32.5 

0.952  ! 

- 

50.20 

35733 

1.68 

—52.5 

— 

25.8 

0.902 

I 

1 

48 

07.00 

31355 

1.92 

— 

-93.3 

40.8 

0.981 

1 

14.80 

29323 

2.05 

-50.2 

— 

43.1 

0.986 

19.50 

28098 

2.14 

-75.2 

25.0 

0.893 

20.30 

27890 

2.16 

-48.3 

— 

26.9 

0.913  1 

29.00 

25623 

2.35 

— 

-66.7 

18.4 

0.783  ! 

33.30 

24502 

2.46 

-43.3 

23.4 

0.874  i 

36.50 

23669 

2.54 

— 

-70.0 

26.7 

0.911  ' 

38.30 

23200 

2.59 

-41.7 

— 

28.3 

0.924 

41.50 

22366 

2.69 

— 

—60.0 

18.3 

0.783  ! 

1 

44.90 

21480 

2.80 

-40.8 

— 

19.2 

0.802 

1 

47.20 

20880 

2.88 

- 

—55.6 

14.8 

0.695 

i 

48.50 

20542 

2.93 

-43.3 

_ 

12.3 

0.608  i 

! 

50.20 

20099 

2.99 

—58.8 

15.5 

0.715  j 

1 

53.90 

19135 

3.14 

-35.9 

— 

22.9 

0.867  1 

1 

55.40 

18744 

3.21 

- 

—63.3 

27.4 

0.917  1 

1 

56.20 

18535 

3.25 

-36.2 

— 

27.1 

0.914 

1 

56.90 

18353 

3.28 

— 

— 57 .5 

21.3 

0.840 

1 

58.80 

17858 

3.37 

-38.2 

— 

19.3 

0.803  1 

! 

49 

00.30 

1 1 4bT 

3.44 

— 

—57.5 

19.3 

0.802  1 

01.10 

17259 

3.49 

-37.2 

— 

20.3 

0.822  j 

01.50 

17154 

3.51 

— 

—52.5 

15.3 

0.710 

04.30 

16425 

3,66 

-37.8 

— 

14.7 

0.690 

' 

05.20 

16190 

3.72 

—55.0 

17.2 

0.756  j 

; 

06.60 

15825 

3.80 

-40.3 

_ 

14.7 

0.690  ' 

1 

07.20 

15669 

3.84 

—54.6 

14.3 

0.678  ! 

08.80 

15252 

3.94 

-38.8 

— 

15.8 

0.725  i 

1 

09.60 

15044 

4.00 

— 

—50.1 

11.3 

0.572 

10.70 

14757 

4.08 

-39.8 

— 

10.3 

0 532 
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Table  XII 

Interference  Pattern  Data.  Mission  15,  Run  8 


■HH 

Range 

ft 

0 

deg 

Max. 

dB 

Min. 

dB 

Diff. 

dB 

“Instantane-  j 
ous  Reflection  i 

Coefficient"  1 

1 

hr 

min 

sec 

17 

08 

26843 

2.24 

-39.3 

12.02 

26034 

2.31 

— 

-64.3 

25.0 

0.893  j 

19.21 

24334 

2.47 

-38.3 

— 

26.0 

0.904 

22.92 

23457 

2.37 

— 

-61.6 

23.3 

0.371  j 

29.25 

21960 

2.74 

-37.8 

— 

23.8 

0.878  ! 

33.84 

20875 

2.32 

— 

-71.8 

0.960  : 

39.53 

19530 

3.08 

-37.0 

— 

34.8 

0.960  1 

42.84 

18747 

3.21 

— 

-51.7 

14.7 

0.690 

45.46 

18128 

3.32 

-35.3 

— 

16.4 

0.737 

48.63 

17378 

3.46 

— 

-64.2 

28.9 

0.930  ; 

16338 

3.68 

-36.3 

— 

27.9 

0.923  ; 

54.79 

15922 

3.78 

— 

-50.8 

14.5 

0.685 

56.88 

15427 

3.90 

-32.7 

— 

18.1 

0.775  1 

58.72 

14992 

4.01 

— 

-49.2 

16.5 

1 

09 

01.81 

14262 

4 22 

-35.4 

— 

13.3 

0.665 

1 

1 

13723 

4.38 

— 

-57.7 

22.3 

0.857 

1 

05.52 

13385 

4.50 

-34.1 

23.6 

0.875 

1 

12727 

4.73 

— 

— 50.7 

16.6 

09.38 

12472 

4.82 

-35.2 

— 

15.5 

10.33 

12200 

4.93 

— 

-46.1 

10.9 

0.556  i 

12.50 

11734 

5.13 

-32.2 

— 

13.9 

13.77 

11434 

5.26 

— 

— 56.9 

24.7 

0.890  ; 

14.58 

11242 

5.35 

-33.1 

— 

23.8 

0.379 

15.35 

11060 

3.44 

-49.2 

16.1 

0.730  1 

16.16 

10869 

5.54 

-31.0 

— 

18.2 

0.780  i 

17.30 

10532 

3.70 

— 

-40.0 

9.0 

0.476 

13.37 

10299 

5.84 

-33.3 

4 7 

4.  1 

0.256  1 

19.79 

10011 

6.01 

— 

-45.7 

10.4 

0.336 

20.23 

9907 

6.07 

-33.0 

— 

12.7 

0.62  1 

21.08 

9706 

— 

-59.4 

26.4 

0.908 

9471 

6.35 

-30.5 

— 

28.9 

0.930 

22.47 

9377 

6.42 

— 

-38.2 

7.7 

0.424 

22.97 

9259 

6.30 

-31.4 

- 

6.8 

0.373 

1 

23.63 

9103 

6.61 

- 

-42.7 

11.3 

0.572 

23.95 

9027 

6.67 

-34.2 

— 

8.5 

0.453 

24.31 

8942 

6.73 

— 

-43.0 

8.8 

0.467 

25.30 

8708 

6.91 

-32.1 

10.9 

0.556 

1 

25.73 

8606 

6.99 

— 

-37.2 

25.1 

0.894 

1 

! 26.13 

3511 

7.07 

-33.3 

— 

23.9 

0.880 

' 26.30 

8424 

7.14 

— 

-43.7 

10.4 

0.536 

1 27.40 

8211 

7.33 

-32.0 

11.7 

0.587 

27.89 

8095 

7.43 

-42.0 

10. 0 

0.520 

7973 

7.55 

-31.0 

— 

11.0 

0.560 

28.31 

7878 

7.64 

- 

-40.1 

9.1 

0.480 

1 

29.33 

7703 

7.81 

-31.4 

8.7 

0.462 

1 j 

7396 

7.92 

— 

—39.5 

8.1 

0.435 

30.27 

7533 

7.99 

-30.7 

— 

8.8 

0.467 

38 


r 


Table  XIII 

Interference  Pattern  Data,  Mission  17.  Run  11 


Time  — GMT 

Range 

ft 

w 

deg 

Max. 

dB 

Min. 

dB 

Diff. 

dB 

“Instantane- 
ous Reflection 
Coefficient” 

hr 

min 

sec 

16 

— 

34 

05.00 

80000 

0.75 

-52.0 

r— — 

41.80 

70500 

0.85 

— 

-83.0 

31.0 

0.943 

35 

13.00 

62750 

0.96 

-50.0 

— 

33.0 

0.955 

49.30 

53500 

1.12 

— 

-76.0 

26.0 

0.904 

36 

19.50 

46300 

1.30 

-46.5 

— 

29.5 

0.934 

44.00 

49500 

1.52 

— 

—71.5 

25.0 

0.892 

59.50 

36000 

1.67 

-44.0 

— 

27.5 

0.918 

37 

16.75 

31800 

1.89 

—66.5 

OO 

0.860 

26.25 

29300 

2.05 

-36.2 

30.3 

0.940 

39.40 

26600 

2.26 

— 

-65.0 

28.8 

0.928 

46.60 

24400 

2.47 

-41.5 

— 

23.5 

0,375 

51.85 

23000 

2.62 

~ 

-70.0 

28.5 

0.925 

57.40 

21600 

2.79 

-33.5 

36.5 

0.970 

38 

01.35 

20600 

2.92 

— 

-60.0 

26.5 

0.909 

09.50 

18500 

3.25 

-37.0 

— 

23.0 

0.867 

21.75 

15750 

3.82 

-56.0 

19.0 

0.800 

24.50 

15000 

4.01 

-36.0 

— 

20.0 

0.816  i 

28.87 

14000 

4.30 

— 

— 59.5 

23.5 

0.875 

30.20 

13500 

4.46 

— 37.7 

_ 

21.8 

0,850  1 

32.00 

13000 

4.63 

— 

-61.0 

23.3 

0.872 

36.25 

12000 

5.01 

—36.5 

— 

24.5 

0.887 

36.92 

11750 

5.12 

-52.0 

15.5 

0.715 

38.30 

11300 

5.23 

-31.0 

— 

21.0 

0.836  j 

39.65 

11200 

5.37 

—53.0 

22.0 

0,852  I 

41.10 

10900 

5.32 

-32.0 

— 

21.0 

0.836  j 

41.90 

10600 

5.68 

— 

-45.0 

13.0 

0.635 

42.60 

10500 

5.73 

-28.5 

16.5 

0.740 

43.50 

10250 

5.87 

- 

—45.8 

17.3 

0.760 

44.90 

9800 

6.14 

-35.8 

— 

10.0 

0.520 

45.40 

9750 

6.17 

— 

-48.0 

12  2 

0.600 

46.75 

9500 

6.33 

-30.0 

— 

18.0 

0.775 

47.35 

9250 

6.50 

-47.0 

17.0 

0.752 

48.30 

8900 

6.76 

-29.5 

— 

17.5 

0.763  i 

1 
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Table  XIV 

interference  Pattern  Data.  Mission  18.  Run  3 


Time  — GMT 


Range 

ft 


03  I 01.80  I 46114 
i 27.00  39240 

I 37.30  36429 

45.30  j 34213 
59.73  ' 30465 


deg 

Max. 

dB 

Min. 

dB 

Diff. 

dB 

1.13 

-53.5 

* 

- 

1.31 

— 

-69.2 

15.7 

1.53 

-48.0 

— 

21.2 

1.65 

— 

—66.5 

18.5 

44.30 

18196 

3.31 

— 

-60.0 

19.5 

47.70 

17244 

3.49 

-39.2 

20.8 

51.20 

16276 

3.70 

— 

-53.7 

14.5 

54.30 

15433 

3.90 

-39.7 

— 

14.0 

35.90 

15000 

4.01 

— 

-62.6 

22.9 

58.30 

14291 

4.21 

-37.0 

— 

25.6 

39.90 

13927 

4.32 



—54.5 

17.5 

02.16 

13224 

4.53 

-36.5 

18.0 

03.40 

13016 

4.62 

— 

-47.0 

10.5 

05.02 

12594 

4.78 

-36.2 

— 

10.8 

06.96 

12090 

4.98 

— 

—66.5 

30.3 

08.35 

11729 

5.13 

-34.5 

32.0 

09.13 

11321 

5.22 

- 

-55.2 

20.7 

09.90 

11326 

5.31 

-36.5 

— 

18.7 

11.40 

10923 

5.51 

-48.7 

12.2 

13.10 

! 10494  I 

5.73 

14.40 

10136  i 

5.92 

‘Tnstantane-  | 
ous  Reflection  j 
Coefficient”  ' 


0.720 
0.840 
0.786 
.823 
.907 
.918 
.865 
.918 
.868 
.865 
.807 
.832 
.685 


Table  XV 

Tabulated  Results  of  All  Runs 


; 

Date 

1974 

Mission 

No. 

Run 

No. 

No.  of 
points 
in  graph 

<7^ -ft. 
predicted 
from  reflection 
coefficient 

a^-ft 

laser 

profil- 

ometer 

error 

\ 

Wind-kts 
A;C  course 
relative 
to  wind 

4/8 

18 

1 

4 

0.60  ± 0.15 

0.63 

1 -8 

UW 

3 

5 

0.65  ± 0.1 

±0.07 

' 0 

U\V 

1 3/8 

17 

9 

2 

= 1.2  (?) 

-0.5* 

+140 

8/ 

10 

2 

-0.6 

-0.5* 

+20 

7.5,LnV44" 

11 

6 

0.45  s 0.05 

-0.45* 

0 

7 DWSO” 

; 

12 

4 

0.50  ± 0.1 

-0.40* 

+25 

6/  i 

3/6 

16 

2 

5 

0.65  ± 0.15 

— 

8/  1 

3 

4 

-1.0 

— 

— 

8/  1 

4 

3 

=*0.6 

— 

— 

3/  1 

i 

5 

4 

0.90  t 0.1 

— 

— 

7/  i 

1 

Avg. 

0.8  ± 0.2 

(0.75 

i 

1 

to  1.0)** 

(±<25) 

1 

3/4 

15 

1 

3 

^0.8 

0.52 

+54 

3 /LAV 14*  i 

j 

2 

2 

*=0.65 

0.55 

-18 

3,U\V26*  i 

1 

5 

5 

0.8  = 0.3 

-O.S 

0 

( 

I 

1 

8 

6 

0.55  t 0.05 

0.64 

-14 

'U\V75*  1 

2/27 

u 

1 

2 

1.08 

14/UW5*  1 

2 

3 

1.18 

14/LAV55* 

3 

2 

1.18 

14,  LAV  ! 

4 

1 

1.18 

13,  DW 

5 

3 

13/LAV27*  ' 

6 

3 

13,  LAV 

Avg. 

0.8  * 0.5 

1.15 

=*-30 

1 

±0.12 

2/23 

13 

0 

2 

*=0.7 

1.58 

—50 

23  DVV  J 

1 

9 

0 

*=1.2 

1.28 

-6 

19,LAV71*  ; 

1 

12 

3 

-o.s 

1.13 

-29 

/xw  ! 

i 

.4vg. 

-0.9  ± 0.3 

1.33 

-32 

1 

•Extrapolated  with  the  same  slope  as  the  wind  speed  (Fig.  16). 
••Visual  estimates  — see  text 


LW  ■ up  wind.  D\V  * down  wind,  X\V  • cross  wind. 
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Table  XVI 

Comparison  of  Remotely  Sensed  and  Measured 
\Vave  Heights  for  the  Qualified*  Runs 


i 

Mission 
i No. 

Run 

Nos. 

firom 

reflection 

ooeff. 

ft. 

Laser 

profil. 

tr 

• error 

t for  runs 
i > 4 pointsl 

Ship 

Tracking 

error 

tall 

daysi 

13 

2.9.  12 

0.9s0.3* 

1.33* 

— 

DDS21 

Optical 

-32 

14 

1 to  6 

0.8i0,3* 

1.13:0.12 

1 

-30 

15 

3 

O.StO.3 

0.3 

0 

0 

" 

S 

0.55t0.05 

0.64 

-14 

’* 

-14 

16 

O 

0.65t0.15 

l0.75 

I — 

** 

— 

■ 

3 

0.9s0.1 

to  1.0i*» 

1 — 

Radar 

— 

IT 

11 

0.45s0.05 

'0.45* 

-0 

*' 

-0 

1 

12 

0,5t0.1 

-0.4' 

-*25 

- -^25 

13 

1 

0.60:0.15 

—3 

DD724 

■' 

-3 

" 

3 

0.63:0.1 

0.65:0.07 

0 

’* 

0 

* Average 

Cl 

rM 

o 

t 

' :o 

1 

> (Exclud* 

' 

1 

1 Mission 

: ■ “Qualified"  means  those  runs  with  4 or  more  points  with  internal  consistency, 

• ■ Average  of  the  runs 

■ From  visual  estimates 

• ■ Extrapolated  with  same  slope  as  wmd  speed. 


Fig.  1 - Continuous  interference  pattern  on  analog  chart  recorder 
Data  of  8 March  1974,  14:31:40  GMT,  Mission  17,  Run  3.  “wvstandard 
deviation  of  water  surface  height'  • 0.65  ft.  Connon  conditions 
for  all  the  1974  data: 

Transmitter-aircraft  altitude  ■ 1000  ft.  unless  othen^ise 
specified.  Shipboard-receiver  height  ■ 50  ft.  Frequencv  ■ 

1.3  GHz.  Polarization  ■ horizontal.  Taken  from  (.2,  Fig.  3''. 
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SIGNAL  LEVEL  - dB 


- 30- 


-40^ 


-50H 


~60H 


-70H 


'80H 


- 90--* 


( I 


111 


/' 


■v 


I 1 1 1 1 1 ! 

4 5 6 7 8 9 10 

ELAPSED  TIME  IN  SECONDS 

(a)  Sharp  miniaium 


Fig.  2 - Different  types  of  interference  mininia,  digital-computer 
plotted.  Data  of  3 March  1974,  Mission  17,  Run  2,  starting  at 
14:32:43  GMT.  - 0.65  ft.  Taken  from  (2,  Fig.  9). 


ELAPSED  TIME  IN  SECONDS 

(b)  Broad  oinimun 


Fig.  2 - Differenc  types  of  interference  minita,  digits I-conputer 
plotted.  Data  of  3 March  1974,  Mission  17,  Run  2,  starting  at 
14:32:43  CMT.  ■ 0.65  ft.  Taken  from  (2,  Fig.  9). 


-30 


-40 


9P  - “IBAST  nVNOIS 


ELAPSED  TIME  IN  SECONDS 
(d)  Small  minima  on  slope  of  a larger  minimum 

Fig.  2 - Different  types  of  Interference  minima,  digital-computer  plotted.  Data  of  8 March 
1974  Mission  17,  Run  2,  starting  at  14:32:48  GMT.  “ 0.65  ft.  Taken  from  (2,  Fig.  9) 


i33iCniM-39NVM 
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Fig.  3 - Comparison  of  radar*  and  optically*tracked  ranges  for  Mission  17  (8  Mar  74). 

The  discrepancy  is  11,000  ft.  USS  JOHNSTON. 


i33J0lW-39NVM 
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1533  1534  1535  1536  1537 

TIME- GMT 

Pig.  4 - Comparison  of  radar-  and  optica lly-tracked  rangea  for  Mission  18  (8  Apr  74).  USS  LAFFEY. 


5 


Pig.  5 - "Instantaneous  reflection  coefficient"  va  grazing  angle.  Mission  13,  Run 
(25  Feb  74).  Manual  optical  tracking.  The  large-circle  points  are  for  the  data 
smoothed  over  1°  grazing-angle  Intervals. 
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54 


Fig.  9 - "Instantaneous  reflection  coefficient*  vs  grazing  angle.  Mission  18,  Run  3 
(8  Apr  74).  Automatic  radar  tracking.  Large-circle  points  are  for  1°  grazing-angle 
Intervals. 
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receiver  on  cliff  at  Chesapeake  Bay.  Transmitter-aircraft  altitude  1000  ft.  Horizontal 
polarization.  Modified  Fig.  7 or  (1). 


INSTANTANEOUS  REFLECTION  COEFFiaENT 


0.7 


Fig.  12  - ’'In«t«nt«n€Ous  reflection  coefficient"  vs  greslng  angle. 

13  Sep  73.  Fixed  receiver  on  cliff  at  Chesapeake  Bay.  Transmitter- 
aircraft  altitude  2000  ft.  Vertical  polarization.  Modified  Fig. 

8 of  (1). 
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REFLECTION  COEFFICIENT 


Fig.  13  - Smoothed  reflection  coefficient  vs  apparent  surface 
roughness,  (o,,3int|»)/X,  compared  to  Ament's  (1953)  theory. 
Mission  15,  Run  8 (4  Mar  74).  Optical  tracking. 

O Plotted  with  " 0.64  ft. 

■ Plotted  with  * 0.55  ft. 

leser  profllometer  • 0.64  ft. 
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REFLECTION  COEFFICIENT 


0 a02  0.04  006  008  0.10  0J2 


Fig.  14  - Smoothed  reflection  coefficient  vs  apperent  surface  rough* 
ness  compared  to  Ament's  theory.  Mission  17,  Run  11  (8  Mar  74') . 
Radar  tracking. 

Data  points  plotted  with  * 0.45  ft. 

Laser  profllometer  a (extrapolated)  « 0.45  ft. 

w 


REFLECTION  COEFFICIENT 


J 


0 0.02  0.04  0.06  0.08  0.10  0.12 

(<r^  SIN  Vf)  / X 


Fig.  15  - Smoothed  reflection  coefficient  vs  apparent  surface  rough- 
ness compared  to  Ament's  theory.  Mission  17,  Run  12  (8  Mar  74). 
Radar  tracking. 

Data  points  plotted  with  o^  - 0.50  ft. 

Laser  profllometer  o^  (extrapolated)  -r  0.40  ft. 

Transmitter-aircraft  altitude  ■ 2500  ft. 


GRAZING  ANGLE -OEG. 


Fig.  17  - Amant'a  theoretical  coherent  field  reflection  coefficient  va 
grating  angle,  with  wave  height  (o^)  aa  a parameter.  Plotted  for  a 
frequency  of  1.3  GHz. 
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.a 


0 I 2 3 4 5 6 7 


GRAZING  ANGLE-DEG 

Fig.  23  - Smoothed  reflection  coefficient  vs  grazing  angle. 
Mission  16,  Run  5 (6  Mar  74).  Radar  tracking. 

Predicted  ■ 0.9  + 0.1  ft. 

Visual  estimates  of  * 0.75  to  1.0  ft. 
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REFLECTION  COEFFICIENT 


C I 2 3 4 5 6 

GRAZING  ANGLE -OEG. 

Fig.  24  - Smoothed  reflection  coefficient  vs  grazing  angle. 
Mission  16,  Runs  3 and  4 (6  Mar  74).  Radar  tracking. 
Predicted  of  points)  ■ 0.8  + 0.2  ft. 

Visual  estimates  of  > 0.75  to  1.0  ft. 


THEORY 

(AMENT) 

ft. 


GRAZING  ANGLE-OEG. 

Fig.  25  - Smoothed  reflection  coefficient  vs  graslng  angle. 
Mission  18,  Runs  1 and  3 (8  Apr  74).  Radar  Tracking. 

Predicted  (Run  1)  ■ 0.60  + 0.15  ft. 

Predicted  (Run  3,  excluding  the  point  at  1.5°)  ■ 0.65  + 0.1  ft 

Laser  profllometer  o„  • 0.65  + 0.07  ft. 


REFLECTION  COEFFICIENT 


Fig.  26  - Smooched  reflection  coefficient  vs  grestng  angle. 
Mission  14,  Runs  1 to  6 (27  Feb  74).  Optical  tracking. 
Predicted  o„(avg.  of  all  points)  = 0.8  + 0.5  ft. 

Laser  profilometer  (avg.  of  all  points)  - 1.15  + 0.12  ft. 
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REFLECTION  COEFFICIENT 


C'j  f-H 
04  cv|  09  00 

00  O ^ ^ O 

CO  «5  O CO  lO 

lO  <0  CO  ^ <0 


i-<  04  lO  CO  rH  01 

^ CO  o>  O UO 

eo  ^ kO  fi  CO 

iO  lO  lO  CO  CO  CO 


QO'^^t^COi09)C^«-^COCOOOO 

'^U06>COcbcOcb^r^iO^^^ 

04CO»HOi*-»CO^iOi-«OICOiOO 

•^^•^^kOiOiOiOCOCOCOCOr^ 


kOO>OOt^UOOOOO 

lO^iMCOvftT-tiOiO^OO 

od6l^'VO^C4CjcOl^ 

OJOJCO^iOiflOOi-^T-* 

COCOCOCOCOCOC^C^C^C* 

0666066666 

C004C0O^C4C0'VrM0l 

c>6^c4d6*owt>6’^ 

f-HC4CO^'^UOOO»-‘i-c 

COCOCOCOCOCOC^C-C^O' 


»-^  lO  f 04  O O 

tH  04  ^ vO  fH 

^ c*  ^ 

UO  O 04  CO 

C«  CO  00  OO  CO  00 


000000 

04  04  ^ ^ 

fM  6 io  ^ 6 

lO  o O 04  04 

C«  00  00  00  CO  00 


Oo^OiOiOO-OOOOi~4 
04  U5  CO  40  04  lO  lO  O C 4 

^<io66^*^6o64HcorM6 

i-t04C0^i0OOi-«04C0'V’^ 

c«t*r*r-r«cooocooooOQoao 

666666666666 

•^OOOOC404OOTf04O 

ii>i^co*Hcidiou9<o^i^6 

O04C0^40i0O*-»04C0C0'«f 

r»c^c*-o-c^c*oooooooooo® 


i-h04C0'^4OC01'®OO»-^04C0 


75 


Appendix  B — Laser  Profilometer  Wave  Heights  (a^ ) 


Mission 

i 

1 

1 

Laser 

Profil. 

1 

Start 

Time 

Stop 

Time 

Laser 

Profil 

Aircraft 

Heading 

No. 

Date 

Run  No. 

GMT 

GMT 

deg. 

ft. 

2/25/74 


2/27/74 


1710:10 

1715:00 

1718:10 

1721:16 

1726:20 

1735:00 

1736:00 

1737:30 

1744:00 

1810:00 

1826:10 

1834:00 

1735:00 

1753:00 

1756:10 

1800:10 

1803:45 

1805:30 

1808:23 

1811:45 

1818:24 

1822:11 

1825:03 


1737:00 

1739:30 

1744:30 

1813:00 

1827:40 

1836:30 

1737:00 

1754:00 

1758:10 

1802:10 

1804:45 

1806:30 

1810:23 

1813:45 

1820:24 

1823:11 

1826:03 


15  3/4/74 


1619:38 

1623:25 

1625:13 

1627:14 

1631:54 

1633:32 

1635:35 

1639:45 

1643:15 

1646:38 


1621:38 

1624:25 

1626:13 

1630:14 

1632:54 

1634:32 

1637:35 

1641:45 

1645:12 

1648:38 
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Appendix  B — Laser  Profilometer  Wave  Heights  (o  ) (Continued) 


Laser 

Profil 

Laser 

Start 

Stop 

•Aircraft 

Mission 

Profil. 

Time 

Time 

Heading 

No. 

Date 

Run  No. 

GMT 

GMT 

deg. 

' = 

12 

f ■ ■ — ■— 

1649:44 

1650:50 

260 

1662:55 

1653-10 

090 

13 

1653:45 

1654-43 

130 

14 

1666:10 

1658:10 

310 

15 

1659:25 

1701:25 

150 

16 

1702:55 

1703:55 

310 

17 

1705:20 

1707:20 

155 

18 

1708:30 

1710:30 

350 

19 

1715:40 

1718:40 

350 

16 

3/6/74 

1 

No  laser  profilometer  data. 

17 

3/8/74 

1 

1422:50 

2 

1427:28 

3 

1435:15 

4 

1438:35 

5 

1445:50 

6 

1451:15 

7 

1459:55 

8 

— 

9 

1512:45 

10 

1521:30 

11 

— 

12 

1530:00 

13 

— 

14 

— 

15 

1544:40 

18 

4/8/74 

1 

2 

3 

4 

Avg. 

APPENDIX  C 


Shipboard  Visual  Observations  of  Sea  State 


Mission  13.  mUlk. 

Estimated  sea  state:  5 

Estimated  wave  height*:  8-10  ft. 

This  was  fairly  consistent  throughout  the  mission  with 
some  decrease  at  end  of.  mission. 

Shipboard  photos:  Available,  but  are  not  tied  in  with  time 
except  for  sequence. 

Aircraft  photos:  Available. 

Mission  14,  2/27/74.  No  observations. 

Mission  15.  3/4/74 

Estimated  sea  state:  0-1.  Winds:  1 to  2 kts. 

Estimated  wave  height:  On  the  order  of  1 ft,  maximum  2 ft. 

Shipboard  photos:  None. 

Aircraft  photos:  Available.  No  whltecaps  on  the  photos. 

Mission  16.  3/6/74. 

Time:  13:04:30  GMT 

Estimated  sea  state:  Smooth,  glassy  seas;  small  localized 
cats  paws. 

Estimated  wave  height:  Local  wind  waves  less  than  1 ft; 
swells  '-I  ft,  occasionally  5 ft. 

Time:  14:06:30  GMT 

Estimated  sea  state:  Small  wavelets  on  waves;  seas  still 
basically  smooth  (on  top  of  the  swell). 

Shipboard  photos:  None. 

Aircraft  photos:  Available. 


* The  observed  wave  height  is  approximately  the  "significant  wave 
height"  ■ “ 4.0o  , where  a Is  the  standard  deviation  of 

height  fluctuations.  ” ” 
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Mission  17.  3/8/74. 


Time  Est . Sea 

State 

Est.  Wave 
Height 

Wind 

Description 

GMT 

ft 

kts 

14:37  2 

2-3 

14 

Occasional  whitecaps 

16:10:30  1-2 

1-2 

7 

Occasional  whitecaps 
sea  state  borderline 
between  1-2. 

17:05 

1-2 

<7 

Only  very  occasional 
whitecaps. 

Estimated  sea  state:  steadily  decreased  during  the  mission. 

Estimated  wave  height:  towards  end  of  mission  the  estimated 
wave  height  was  still  1-2  ft.  The  wind  dropped  steadily 
during  the  mission.  At  the  end  of  the  mission,  the  waves 
were  still  decreasing  and  probably  closer  to  1 ft  than  2 ft. 
By  evening  waves  were  less  than  1 ft . 

Shipboard  photos:  Available. 

Aircraft  photos:  Available. 

Mission  18,  4/8/74. 


Time 

Est . Wave 

Description 

GMT 

Height 

ft 

11:25 

1 

Only  occasional  whitecaps. 

11:56 

perhaps 

1 

No  white  caps. 

12:39 

1 to  1 1/2 

Only  occasional  whitecaps. 

13:41 

1 to  1 1/2 

Only  occasional  whitecaps. 

Shipboard  photos:  Available. 

Aircraft 

photos:  Available. 
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APPENDIX  D 


Shipboard  Wind  Observations 
Anemometer  Is  approximately  113  ft  above  water. 
Mission  13.  2/25/74. 


wind  Data 


Mission  17.  3/8/74. 


Time 

Speed 

Wind 

GMT 

Direction 

kts 

true  - deg. 

1355 

15 

286® 

1 

1420 

17 

290® 

I 

1433 

(16)* 

i 

1453 

14 

300® 

I 

1500 

12 

300® 

1507 

13 

300* 

1516 

10 

315® 

1521 

11 

314® 

1541 

12 

313® 

1557 

10 

309* 

1603 

11 

330* 

1604 

11 

314® 

1608 

(9)* 

1614 

8 

313® 

1621 

8 

310® 

1623 

7 

318® 

1629 

7 

321® 

1635 

7 

320® 

1638 

- 

320® 

1650 

(6)* 

1656 

6 

316® 

Mission  18.  U/8/7^. 

1210  Orr  22  kts  039® 


il 


* ■ Interpolated. 


APPENDLX  E 


Forward-scatter  remote  sensing 
15  April  1978 

From  further  analysis,  it  now  appears  that  forward  scattering  may  be 
capable  of  determining  several  characteristics  of  ocean  waves. 

1.  Wave  height 

Below  an  apparent  roughness  of  0.1,  Ament's  (1953)  theoretical 
curve  is  recast  as  a family  of  curves  vs  grazing  angle  (0),  with  wave 
height  (a  ) as  a parameter.  Plotting  the  measured  reflection  coeffi- 
cients of  the  L-band  data  on  this  graph  then  yields  wave  height  directly. 

2.  Surface  correlation  length  along  line-of-sight 

Above  an  apparent  roughness  of  0.1,  two  plane  wave  theories,  one 
of  acoustical  scattering  from  a surface  with  exponential  height  statis- 
tics [Boyd  and  Deavenport,  1973],  and  the  second  of  EM  scattering  from 
a particular  non-Gaussian  surface  [Brown  and  Miller,  1974],  both  follow 
the  over-ocean  experimental  data  [Beard,  1961].  Boyd  and  Deavjnport 
also  pointed  out  the  agreement.  Assuming  that  the  agreement  is  meaning- 
ful, even  with  non-Gaussian  surfaces,  led  to  the  realization  that 
shadowing  of  waves  by  other  waves  is  a mechanism  by  which  a non-Gaussian 
surface  can  be  created  from  a Gaussian  ocean.  The  probability  distri- 
bution of  the  non-shadowed  (i.e.,  illuminated)  portion  of  the  ocean 
surface  has  been  derived  by  Wagner  [1967],  and  is  non-Gaussian.  In 
addition,  the  effective  wave  height  of  the  non-shadowed  surface  is  less 
than  the  physical  wave  height  of  the  ocean  surface.  Both  these  factors 
would  raise  Ament's  curve  upward  towards  the  ocean  data  "experimental 
curve. " 


Wagner’s  shadowing  theory  is  a function  of  a , tanv,  and  L , 
where  L is  the  surface  correlation  length  along  the  direction  of 
the  propagation  path  (from  transmitter  to  receiver) . The  unknown  a 
can  be  found  as  above,  leaving  only  as  the  one  unknown  to  be  solved 

for.  Wagner's  expression  needs  to  be  tested  for  its  region  of  appli- 
cability where  values  of  L produce  best  agreement  with  experimental 
data.  ^ 


3.  Surface  correlation  length  transverse  to  line-of-sight 

The  fact  that  Beckmann's  [1967]  spherical  wave  theory,  evaluated 
for  wave  tank  conditions  [Beard,  1967],  follows  the  over-ocean  data  also 
seemed  remarkable.  Calculations  made  of  the  sphericity  parameters  of 
the  wave  tank  conditions  indicate  that  the  wave  height  (in  the  vertical 
plane)  is  insufficient  to  produce  the  observed  sphericity  effects. 
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However,  if  one  looks  down  on  the  ocean  at  scattering  in  azimuth,  or 
in  the  horizontal  plane,  the  ocean  correlation  length  transverse  to 
line*^f “Sight  path  CLy)  is  sufficient  to  produce  the  observed 
sphericity  effects,  CLy  can  be  thought  of  approximately  as  the  wave™ 
crest  lengths.)  ^ 

Beckmann's  theory  is  an  explicit  function  of  a , 0,  L , R,  and 
where  R;  is  the  range,  or  distance  between  trnasmitter  and  receiver. 
It  also  seems  to  be  an  implicit  function  of  L , and  thus  has  the  capa- 
bility of  determining  L^.  ^ 

4,  Proposed  remote  sensing  experiment 

A millimeter  wave  experiment  at  David  Taylor  Model  Basin  is 
proposed  for  FY-79  to  exploit  the  different  functional  dependencies  of 
these  theories  to  determine  the  wave  characteristics  a , L , and  L . 
The  region  of  applicability  of  each  theory  has  first  to  be’^determiXed . 

All  of  these  measurements  for  different  functional  dependence 
must  be  made  under  controlled  conditions  and  with  adjustable  water 
wave  directions  (to  vary  ^ and  with  constant  o^,  etc.).  Conse- 
quently, a wave  basin,  as  at  David  Taylor  Model  Basin,  is  required. 

The  two-dimensional  surface  correlation  function  must  be  measured. 

In  order  to  achieve  apparent  surface  roughnesses  between  0.1 
and  0.3,  where  the  various  theories  differ,  it  is  necessary  to  use  a 
short  wavelength,  3.3  mm  (90  GHz). 


C.  I.  Beard 


